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ABSTRACT 



A pipeline technique by Loomis and Sinha has been applied to the design of recur- 
sive digital filters. Recursive digital filters operating at hitherto impossibly high rates can 
be designed by this technique. 

An alternate technique by R. Gnanasekaran allows high speed implementation using 
the state- space structure directly. High throughput is also achieved by use of pipelined 
multiply-add modules. The actual hardware complexity will depend upon the number 
of pipeline stages. 

These techniques are used for the design of the HR notch filter and finally, a com- 
parison of the performance and complexity of these two techniques is presented. 



C.l 

TABLE OF CONTENTS 

I. INTRODUCTION 

A. DIRECT FORM STRUCTURE REALIZATION 

B. STATE SPACE REALIZATION 

C. OVERVIEW OF THESIS 

II. HR NOTCH FILTER DESIGN 

HI. REALIZATION OF HR NOTCH FILTER BY USING THE 
LOO.MIS SINHA METHOD 

A. PIPELINE RECURSIVE DESIGN 

B. SECOND ORDER HR NOTCH FILTER: P = 0 

C. AUGMENTED HR NOTCH FILTER: P = I 

D. AUGMENTED HR NOTCH FILTER: P = 2 

E. AUGMENTED HR NOTCH FILTER: P = 3 

F. STABILITY TESTING 

1. Condition for Stability of the Original HR Notch Filter: p = 0 

2. Condition for Stability of the Augmented HR Notch Filter: p = 1 

3. Condition for Stability of the Augmented HR Notch Filter: p = 2 . . . . 

4. Condition for Stability of the Augmented HR Notch Filter: p = 3 . . . . 

G. FREQUENCY RESPONSE 

1. Second Order HR Notch Filter: p = 0 

2. The .Augmented HR Notch Filter: p = 1 

3. The Augmented HR Notch Filter: p=2 

4. The Augmented HR Notch Filter: p=3 

H. UNIT SAMPLE RESPONSE 

1. Second Order HR Notch Filter: p=0 

2. The Augmented HR Notch Filter: p = 1 

3. Ihe Augmented HR Notch Filter: p=2 

4. The Augmented HR Notch Filter: p=3 

I. PIPELINE IMPLEMENTATION 

1. Second Order HR Notch Filter: p = 0 



. 1 

. 1 

. I 

. 1 

. 3 

. 8 

. 8 

11 

12 

15 

17 

20 

20 

21 

'>•> 

23 

26 

26 

26 

26 

26 

26 

26 

26 

26 

34 

34 

34 



IV 



2. The Augmented HR Notch Filter: p = 1 34 

3. The Augmented HR Notch Filter: p = 2 36 

4. The Augmented HR Notch Filter: p = 3 36 

J. SUMMARY 39 

IV. GNANASEKAR^AN’S STATE SPACE REALIZATIONS 41 

A. STATE SPACE REPRESENTATION OF A NTH ORDER RECURSIVE 

DIGITAL FILTER 41 

B. WRITING THE STATE SPACE EQUATIONS 42 

C. MODIFIED STATE SPACE EQUATION 43 

D. SECOND ORDER HR FILTER STATE SPACE REALIZATIONS 44 

1. L = 1 45 

2. L = 2 46 

3. L = 3 49 

4. L = 4 52 

E. PIPELINED IMPLEMENTATION 56 

1. L = 1 56 

2. L = 2 58 

3. L = 3 • 58 

4. L = 4 62 

5. Example for L = 1: 62 

a. Ka = 0. Km = 1 62 

6. Example for L = 2: 62 

a. Ka = 0. Km = 2 64 

b. Ka = 1. Km = 0 64 

7. Example for L = 3: 67 

a. Ka = 0. Km =3 67 

b. Ka = 1. Km = 1 67 

8. Example for L = 4: 71 

a. Ka = 0, Km = 4 71 

b. Ka = 1, Km = 2 71 

c. Ka = 2. Km = 0 71 

F. SUMMARY 74 



V 



V. CONCLUSIONS 76 

A. SUMMARY 76 

B. COMPARISON 76 

APPENDIX A. PROGFLAM AND DATA OF FREQUENCY RESPONSE 78 

A. PROGIUAM FOR FREQUENCY RESPONSE CALCULATION 78 

B. MAGNITUDE, PHASE OF SECOND ORDER HR NOTCH FILTER: 

P = 0 88 

C. MAGNITUDE, PHASE OF THE AUGMENTED HR NOTCH FILTER: 

P=1 91 

D. MAGNITUDE, PHASE OF THE AUGMENTED HR NOTCH FILTER: 

P = 2 93 

E. MAGNITUDE, PHASE OF THE AUGMENTED HR NOTCH FILTER: 

P=3 96 

APPENDIX B. PROGIGAM AND DATA OF DIFFERENCE EQUATION ... 99 

A. PROGR.AM FOR DIFFERENCE EQUATION CALCULATION 99 

B. INPUT, OUTPUT DATA OF SECOND ORDER HR NOTCH FILTER: 

P = 0 108 

C. INPUT. OUTPUT D.ATA OF THE AUGMENTED HR NOTCH FIL- 
TER: P= 1 116 

D. INPUT. OUTPUT DATA OF THE .AUGMENTED HR NOTCH FIL- 
TER: P = 2 124 

E. INPUT. OUTPUT DATA OF THE AUGMENTED HR NOTCH FIL- 
TER; P= 3 132 

LIST OF REFERENCES 141 

INITIAL DISTRIBUTION LIST 142 



VI 



LIST OF TABLES 



Table 1. b 's for/? = 0.1 5 10 

Table 2. a s for/> = 0,1,. ..,5 11 

Table 3. a 's for;? = 0,1,..., 3 11 

Table 4. LOOMIS, SINHA PIPELINE RECURSIVE DESIGN 39 

Table 5. STATE SPACE REALIZATION SU.M.MARY 74 

Table 6. SUMMARY OF REALIZATIONS 77 



LIST OF FIGURES 



Figure 1. Poles and Zeros Location 6 

Figure 2. Second Order HR Notch Filter; p=0 7 

Figure 3. Second Order HR Notch Filter: p = 0 13 

Figure 4. The Augmented HR Notch Filter: p = 1 14 

Figure 5. The Augmented HR Notch Filter: p =2 16 

Figure 6. The Augmented HR Notch Filter; p=3 19 

Figure 7. Condition for Stability of the Original Filter 21 

Figure 8. Condition for Stability of the Augmented Filter; p = 1 23 

Figure 9. Condition for Stability of the Augmented Filter: p = 2 24 

Figure 10. Condition for Stability of the Augmented Filter: p =3 25 

Figure 11. .Magnitude Response of Second Order HR Notch Filter; p = 0 27 

Figure 12. .Magnitude Response of the Augmented HR Notch Filter: p = 1 .... 28 

Figure 13. .Magnitude Response of the .Augmented HR Notch Filter: p = 2 .... 29 

Figure 14. .Magnitude Response of the Augmented HR Notch Filter: p = 3 .... 30 

Figure 15. Output Sequence of Second Order HR Notch Filter: p = 0 31 

Figure 16. Output Sequence of The Augmented HR Notch Filter: p = 1 32 

Figure 17. Output Sequence of The Augmented HR Notch Filter: p = 2 33 

Figure 18. Output Sequence of The Augmented HR Notch Filter: p = 3 35 

Figure 19. Second Order HR Notch Filter; p = 0 with Ka = 1, Km=0 36 

Figure 20. The Augmented HR Notch Filter: p = 1 with Ka = 1, Km= 1 37 

Figure 21. The Augmented HR Notch Filter: p = 2 with Ka= 1, Km=2 38 

Figure 22. The Augmented HR Notch Filter: p = 3 with Ka= 1, Km=3 40 

Figure 23. Second Order HR Filter; L = 1 47 

Figure 24. The Augmented HR Filter: L = 2 50 

Figure 25. The Augmented HR Filter: L = 3 53 

Figure 26. The Augmented HR Filter: L = 4 57 

Figure 27. Heart of the Pipeline Recursive Filter Realization 58 

Figure 28. Second Order HR Filter; L = 1 with Ka, Km 59 

Figure 29. The Augmented HR Filter: L = 2 with Ka, Km 60 

Figure 30. The Augmented HR Filter; L = 3 with Ka, Km 61 

Figure 31. The Augmented HR Filter; L = 4 with Ka, Km 63 

viii 



Figure 32. Heart of More General Pipeline Filter Realization 64 

Figure 33. Second Order HR Filter: L = 1 with Ka = 0, Km= 1 65 

Figure 34. The Augmented HR Filter: L = 2 with Ka = 0. Km=2 66 

Figure 35. The Augmented HR Filter: L = 2 with Ka= 1, Km=0 68 

Figure 36. The Augmented HR Filter: L = 3 with Ka = 0, Km= 3 69 

Figure 37. The Augmented HR Filter: L = 3 with Ka= 1, Km= 1 70 

Figure 38. The Augmented HR Filter: L = 4 with Ka = 0, Km = 4 72 

Figure 39. The Augmented HR Filter: L = 4 with Ka= 1, Km=2 73 

Figure 40. The Augmented HR Filter: L = 4 with Ka=2, Km=0 75 



IX 



I. INTRODUCTION 



A. DIRECT FORM STRUCTURE REALIZATION 

A scheme for high speed recursi\'e digital filter realization was suggested by H. H. 
Loomis and B. Sinha [Ref 1]. The main aspect of the design is in modifying the filter 
difference equation so that the recursive part requires only the Lth past output instead 
of the immediate past output. This removes the tight recursive requirement thus enabling 
pipelining of multipliers and adders in the direct form implementation resulting in high 
throughput. One of the difficulties with this method is that the modification introduces 
extra poles (not necessarily stable) and zeros in the filter transfer function which cancel 
each other in the infinite arithmetic implementation. In practice due to the finite word 
length efiect this cancellation may not occur. This may result in an unstable realization 
even if the original filter were a stable one. However, if L is chosen to be large enough, 
then it has been proven that the new poles introduced will also be within the unit circle 
thus guaranteeing stability even if the pole-zero cancellation does not occur. This 
method is called the direct form structure, or the scalar structure. 

B. STATE SPACE REALIZATION 

Another method by R. Gnanasekaran is an alternate technique which relaxes the 
recursive nature of the state equation (state space realization) [Ref 2|. That is, the 
present state vector computation does not require the immediate past state vector. The 
proposed technique allows high speed implementation using the state-space structure 
directly. High throughput is achieved by use of pipelined multiply -add modules. The 
actual hardware complexity will depend upon the method of pipelining, though very high 
speed is possible by inc’easing the number of states in the pipelining. This method im- 
plements A^ instead of A. the system matrix. The statespace realization will be stable 
if the original filter is stable. Furthermore, the A^ filter has improved tolerance to coef- 
ficient error. 

C. OVER^TE^y OF THESIS 

Chapter II describes the HR notch filter. The transfer function and the coefficients 
that are appropriate for the HR notch filter are developed. 

Chapter III presents the realization of the high speed HR notch filter by using the 
pipeline technique of Loomis and Sinha. Stability testing, frequency response and unit 



1 



sample response are described. System designs for pipelined multiply-add modules are 
included. 

Chapter IV presents Gnanasekaran's state space realization which includes a modi- 
fied state space equation. System designs for pipelined multiply-add modules are de- 
scribed. 

Chapter V contains a comparison of the results of the Loomis/Sinha pipeline re- 
cursive design and the Gnanasekaran state space realization design. 

Appendix A contains the program and data for frequency response analysis of the 
HR notch filters in Chapter III. 

Appendix B contains the program and data for unit sample response analysis of the 
HR notch filters in Chapter HI. 
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II. HR NOTCH FILTER DESIGN 



From Rowe, [Ref. 3: p. 72] the HR notch filter has the transfer function 



H{2) = 



6qZ + b]Z -f ^0 
Oq2^ + b^z + «2 



( 2 . 1 ) 



Dividing the numerator and denominator of (2.1) by Oq, we have 



H{z) = 



^0 2 , ^1 , ^0 
77~ z + “77“ 2 4- “77“ 
^ ^ 

2 ^1 

' + — ^ ■’ H ^ 

- ^ flo - ^ ao 



Let 



«0 



„ 4>, „ , a, 

5o, ^ = 5,. and 



Ag. Therefore the transfer function is 



m = 



BrfZ^ + 5jC + Bq 

z- + 5,r -r^o 



( 2 . 2 ) 



(2.3) 



From (2.1) and (2.3) we can see that 

1. 5, in the numerator, must be equal to 5, in the denoininator; and, 

2. the first Bg in the numerator must be equal to the last Bg in the numerator. 

Rewriting (2.3) as 



2 ^1 

H{z) = r 2 , (2.4) 

(j + 5)Z T Aq) 



the transfer function H{z) is of the form 



Boiz-n,)(z- 
{z-d,){z-d.) 



(2.5) 



Let the zeros be placed at Ic-®, and let the poles be placed at Therefore 

the transfer function is 



3 



(2.6) 



B^iz - \e-^^){z - 



(r - e~n{z - v''^o 



le-^"z+l] 



-je. 






[2^ - JAq - V^o + JA ] 



H{z) = 



2 

B^\?-l{^^ ^ - - )z+l] 



[2 - 2V^o ( 






)z + Aq] 



Hiz) = 



Bq\_z — 2 cos 6z +1J 



[r — cos + yio] 



5qz —2BqCOS$z + Bq 

n{z) = -1 — = 

2 — 2^’Aq cos 9z + Aq 



According to (2.3) we can see that 
—2Bq cos $z —- —2JAq cos 6z 



(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

(2.11) 



which means 

Bo = ^- (2.12) 



For our e.xample of the HR notch filter we want the zeros to be placed at 

j^±/0.4467110339z- ^ ^^±/e 



(2.13) 



and we want the poles to be placed at 

0_99^±/0.4467U0339. ^ 



(2.14) 



We can see that 



0 = 0.4467110339;: (2.15) 

^0 = 0.9801 (2.16) 



Bq = JAq =0.99 



(2.17) 
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This realization is shown in Figure 1 on page 6. 
From equation (2.1 ) 



r,, ^ B^z- -2Bo cos dz + Bq 

Ji(z) = 

r — 2^;Aq cos $z + Aq 



Substituting 9, Bq and A^. we then have 

0.99r^ - 2(0.99) cos(0.44671 lQ339re)z + 0.99 
2 ^ - 2(0.99) cos(0.4467110539;r)z + 0.9801 



H{z) = 



0.99---0.32993Z + 0.99 
z*- 0.32993r-f0.9S01 



//(z-') = 



0.99 - 0.32993z~' 0.99z~^ 

1 -0.32993z"' +0.9801r~^ 



Because we have the transfer function [Ref. 1: p. 273] 



= -TU,r 



\ - b,z~^ - b,z~^ 



from (2.21) and (2.22) we can see that all of coefficients are 



«o = 0.99 
= -0.32993 
= 0.99 
Z), =0.32993 
b. = -0.9801 . 






>'(--) 

A'(z) 



>'U) _ 0.99 - 0.32993z~* + 0.99z~^ 

1 - 0.32993z"‘ + 0.9801z"- 

F(z)[l - 0.32993z"’ + 0.9801z'^] = A'(z)[0.99 - 0.32993z"’ + 0.99z”^] 



(2.18) 

(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 

(2.24) 

(2.25) 

(2.26) 



F(z)-0.32993z"'l'(r)4-0.9S01z‘-y(z) = 0.99A'(z)-0.32993z 'A'(r) + 0.99z -A'(z) (2.27) 




Figure 1. Poles and Zeros Location 

y(r) = 0.99A'(r) - 0.32993r"’.V(2) + 0.99z"^X(z) + 0.32993z"’ F(z) 

7 (2.28) 
-0.9S01Z ^F(‘) 

y{n) = 0.99.r(«) - 0.32993.v(/j - 1) + 0.99;r(rt - 2) + 0.32993>-(« - 1) 

-0.9S01j(rt- 2). 

This is the basic second-order recursive filter (p = 0). The number of p will be ex- 
plained later in Chapter III and the realization is shown in Figure 2 on page 7. 

We can design this HR notch filter to be a high speed HR notch filter by using 
pipeline techniques so that it can be operated at hithero impossibly high rates. In 
Chapter HI we will apply Loomis 'Sinha method and in Chapter IV Gnanasekaran's 
state space realization. A comparison of the results of these two techniques will be 
shown in Chapter V. 
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~~u — = 




Unit Delay 
Pipeline Multiplier 



-e- = 

x(n) 

y(n) 



Pipeline Add Unit 

Input 

Output 



Figure 2. Second Order HR Notch Filter: p = 0 
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III. REALIZATION OF HR NOTCH FILTER BY USING THE 
LOOMIS/SINHA METHOD 



A. PIPELINE RECURSIVE DESIGN 



From Chapter II we have the second order HR notch filter that is described by the 
transfer function 



^(2 ) = n ZT 

(I - r>,z ’ - b 2 Z ) 



(3.1) 



and the difference equation 



y{n) = b^y{n — I) + b^yin — 2) + aQx{n) + a,jr(n — 1) + 02 x{n — 2). 



(3.2) 



The computation of r(/?) requires the immediate past output — I) . This filter can 
be developed into high speed HR notch filter by using the pipeline technique of Loomis 
and Sinha [Ref 1]. The main aspect of the design is in modifying this filter difference 
equation so that the recursive part requires only from the Lth past output instead of the 
inxmediate past output. This means the speed limits can be overcome by modifying this 
difference equation so thatj-(/j) depends not on the immediate past output, — 1), but 
on — ^) for some L ' 1. This can be accomplished as follows: 

For L = 1 the difference equation is 



We can see that this equation is a 2-order difference equation. 

From [Ref 1: p. 274] the difference equation has (p -f- m)^^'-order where m is the order 
of the difference equation and p is the pipeline delay used to overcome the speed limit. 
Therefore equation (3.3) has the number of p = 0. We can call this equation the second 
order HR notch filter with p = 0. 

Delaying equation (3.3) by 1 yields 

y{n — 1) = biy’{n — 2) 4- — 3) -f- aQx{n — 1) -f a^x{n — 2) -t- a 2 x{n — 3). (3.4) 

For L = 2. putting (3.4) into (3.3) results in the difference equation 



;•(«) = biy{n — 1) + b2}in — 2) -h Oov(/i) + + a^x{n -1)4- a2x(n — 2). 



(3.3) 
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(3.5) 



y{n) = — 2) + b 2 }i>i — 3) + ao.x(« — 1) + a^x{n — 2) + a 2 x{n — 3)] 

+ byy{H — 2) + (3o.v(rt) + <7].x(/2 — 1) + — 2) 

y{n) = blv{n - 2) + b^b^yin - 3) + ciQb^x{n - 1) + a^b^x{n - 2) + a 2 b^x{n - 3) 
+ b2y{n — 2) + aQx{n) + a,x(« — 1) + a2x{n — 2) 

>-(m) = (b^ + b2)y{n - 2) + b^b2y{n - 3) + aox(«) + (a^b^ + a,)x(« - 1) 

+ ia^b■^ + a2)x{n — 2) + a2bix{n — 3). 



We see that the computation of>'(«) requires only the second past output y(rt — 2) 
instead of immediate past output y(«— 1), and we can use this technique for the 
higher-order of L . 

We see that equation (3.7) is a 3-order difference equation and has p = 1. Thus we 
can call this filter the augmented IIR notch filter for p= 1. 

The z transform of equation (3.7) is 

} = (/fif - 1 - ^>2); ^ } -f- £>,^2- ^ T-l- -I- (ugf’] -I- a,)z ^ 

-f- (U]^i " 1 - ^2)2 ^X -l- (12b 

}' ^0 + (up^i -r aQz ’ -I- {a^b^ -f 02)2 ^ -I- ^2^-^!- ^ 

I - (b^, + b2)z-^ - b,b2Z-^ ^ 



17,/ -u (^0 + '+«2‘ ^)(I + ') 

^ ) = n 3 rr~ • 

(1 - b^z - b2Z ■)(! + b^z ) 



(3.10) 



That means the transfer function of the augmented IIR notch filter, p= 1. is equal 
to the transfer function of the second order IIR notch filter, p = 0, multiplied and divided 
by the same polynomials. 

We can design the higher-order difference equation by starting from a second-order 
difference equation. That means we can design the higher order, p-augmented HR notch 
filter by starting from our second order HR notch filter, p=0, as follows; 

From Chapter II the second order HR notch filter transfer function is the original 
transfer function 



«(2-') 



(gp -b ( 7 ,z ‘ -I- G 2 Z ^) 

(1-V’"'-V') ‘ 



(3.11) 
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Multiply the numerator and denominator of the original transfer function by the 
same polynomial equation. This can be thought of as equation (3.12) As developed in 
[Ref 1; p. 2S1], the transfer function of the p-augmented filter is as shown in equation 
(3.13) We know from other considerations that the transfer function of the higher order 
realization must in fact be the same as the original and therefore must have the same 
poles and zeroes in the z - plane. 

The coefficients of the a polynomial depend on the original and b^, and p is de- 
termined by the number of stages in the pipeline multiply and add units. 



(uq + ^)(I + ’ + ••• + 

(1 - - b 2 Z~\\ + afz"' + - + af z"0 



(3.12) 



H'(z-') 






, ,(P) -P-2 



/,(P) ,-P-2 

- - Opj^2^ 



(3.13) 



The values of bf, af\ and a!f are given in Tables 1, 2 and 3 [Ref 1: pp. 278-279]. 



Table 1. b 's for p = 0.1.. ...5 





Coefficient of Z) 


Coefficient of 


p 




b^l. 


0 


bs 




I 


h] + bj 


bA 


2 


+ 2bA 


bA + bi 




+ 3bjb2 bj 


bA + 2bA 


4 


b\ + + 7>b^b\ 


b\h2 + l>b\b\ + b\ 


5 


b\ + 5b\b. + 6b\b\ + b\ 


b]b, -I- 4b'A + 3i,i>| 
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Table 2 . a s for p = 0.1.. ...5 



p 






o.f 


ly.p 






0 


1 


0 


0 


0 


0 


0 


1 


1 


b, 


0 


0 


0 


0 


2 


1 


b, 


b\ + bi 


0 


0 


0 


5 


1 


bx 


b\ + b. 


b] + 2b 2 


0 


0 


4 


1 


bx 


bx + ^2 


b\ + 2b^b2 


b\ + l>b\b2 + b\ 


0 


5 


1 


bx 


bx + ^2 


b\ + 2b ^bi 


b\ + WA + b\ 


+ AbA + IbM 



Table 3. a 's for p = 0.1 3 



P 


Ci'f) 

'’*0 


a'f 


aV 


flf 


a'p 




0 


«o 






0 


0 


0 


1 


«0 


c7pi>j -r 0.^ 


rt,Z), 4- fl;. 


a.b^ 


0 


0 




^0 


*r c7j 


aj^b\ + b.) 
+ £7]i*l 4" flj 


a,{b\ + ^2) 

-j- 


a,(^f 4- b^ 


0 






aA + «i 


axkb'xA- />;) 
~h Clih\ ~l“ CI2 


+ (ixib\ + b^ 
4- a2b^ 


-f ci^{b\ -f 'Ib^b'^ 

-f ^2(^1 + ^2) 


4- a^Ax + 2^|/\) 



B. SECOND ORDER HR NOTCH FILTER: P = 0 

From equation (3.13) in section A, we let p = 0. This equation becomes: 



-K («0 + «l- +«2-" ) 

bn^ ) - 


(3.14) 


(1 - yjz - b2Z ") 




N 

II 


(3.15) 


7(2) (flo + «lZ ’+^22 


(3.16) 


(1 — ^i2~’ — ^2^”*) 


.Manipulating (3.16) will result in the difference equation, (3.20). 




}'(z)[l - b^ 2 ~^ - b^z'^l = A'(2)[flo + OlZ”' + « 22 ”‘] 


(3.17) 
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F(z) - ^, 2 -' Y{z) - hz-h\z) = a^X{z) + a,z-'X{z) + a^z'^ATz) (3.18) 

Y{z) = aoA'(-) + a,z-\\\z) + a^z-^X'iz) + b,z~^ }'(z) + b^z-^Xz) (3.19) 

>•(/?) = %r(A?) + a,.v(« - 1) + « 2 -'^(» - 2) + b^y{n - 1) + b^v{n - 2). (3.20) 

This realization is shown in Figure 3 on page 13. 

From chapter II we know that for our example notch filter, that the coefficient 
values are 



no = 0.99 
ai = -0.32993 
«2 = 0.99 
= 0.32993 
= -0.9801 . 



The transfer function is 



//'(z-') = 



(c7q -{- f 

(l-Z),z-'-Zv"1 ' 



The difference equation is 

v(«) = no.v(«) + n,.v(n - 1) + n 2 -^(/' ~ 2) + b^y{n - 1) + b^-{n - 2). 



(3.21) 



(3.22) 



(3.23) 



C AUGMENTED HR NOTCH FILTER: P = 1 

From equations (3.12) and (3.13) in section A, we let p = 1. Those equations be- 
come; 



//'(z-’) = 



(Cq "f * "f (^ 2 ^ "f *) 
(1 - Z?,z"' - b 2 Z~^)(l -1- a^i'^z"') 



//'(z-’) = 



, 0 ) 



, ( 1 ) -1 , ( 1 ) -2 , ( 1 ) 
+ a\'z + 2 + a] z 

i r (1) -2 ,Ti) -3 

I — O 2 2 — O 2 2 



-3 



(3.24) 



(3.25) 



/f'(z-‘) 



iM 

A'(z) 



(3.26) 
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X 




Figure 3. Second Order HR Notch Filter: p = 0 






X{z) 



1 - b 






Manipulating (3.27) will result in the difTercnce equation, (3.31). 

}-(r)[l - = X(z)[c4‘^ + 

y(r) - b^'^-^Xz) - 6]'>z-^ Viz) = 4'lV(z) + af>z-'A'(z) + 4^^z~^X{z) 

+ a^Pz~^X{z) 

Y{z) = 4'^'V(z) + a\%-\\Xz) + «l'^ 2 "lV(z) + 4^^z-^X{z) + F( 2 ) 

+ 6<’>z-'}T2) 

^•(rt) = c^^^x{n) + a{'^.r(rt — 1) + a2^x{n — 2) + a3'^jr(n — 3) + b2^y{n — 2) 
+ - 3). 



(3.27) 

(3.28) 

(3.29) 

(3.30) 

(3.31) 



This realization is shown in Figure 4 on page 14. 

From Tables 1, 2 and 3 [Ref. 1: pp. 278-279], we have 
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X 




Figure 4. The Augmented HR Notch Filter: p = 1 



a</> = ^, =0.32993 
4’> = no = 0.99 
a|'^ = no^i + a, = -0.0032993 
= + fl2 = 0.8811462 

a<'> = a2^>i =0.3266307 
b^^'> = b\ + b2 = -0.^inA(>2 
= = -0.32336439. 



The transfer function is 



//'(2-‘) = 



,0) 






(l),-3 



+ a\'z + a\‘’z + a\'z 



1 _ ^,< 0,-2 _ 



The difference equation is 



(3.32) 



(3.33) 
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(3.34) 



y{n) = + a\'^x{u - 1) + - 2) + 4'^x{n - 3) + - 2) 

+ ^<'>v(n-3). 

D. AUGMENTED HR NOTCH FILTER: P = 2 

From equations (3.12) and (3.13) in section A, we let p = 2. Those equations be- 
come: 



1 _ (gp + a\Z * 4- a2Z ^)(1 + ‘ + a.f'^z ^) 

(1 — — b2Z~~){\ + -I- a2^^z”^) 



(3.35) 




(3.36) 




(3.37) 




(3.38) 



Manipulatiiig (3.38) will result in the diflerence equation, (3.42). 





(3.40) 




(3.41) 



y{n) = a^^x{n) -I- af^x{n — 1) 4- a^^x{n — 2) + 4^^x(n — 3) -f a^^x{n — 4) 
+ ^fV(«-3) + ^f>-(«-4). 



(3.42) 



This realization is shown in Figure 5 on page 16. 

From Tables 1, 2 and 3 [Ref. 1: pp. 278-279], we have 
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Figure 5. The Augmented IIR Notch Filter: p = 2 



af> = Z), = 0.32993 
a?> = Z>f + Z?2 = -0-S7 12462 
4 ') = «o = 0-99 
af^ = oo^i + fli = -0.0032993 
02^ — ^(^1 "H ^ 2 ) ^1^1 "1“ ^2 ~ 0.01861246 

af = + b2) + ^2^1 = 0.61408096 

a<-> = a2(^f + ^ 2 ) = -0.86253373 
+ 2Z>,^2 = -0.61081465 
) = Z9-^2 + ^2 = 0.8539084 . 



The transfer function is 



, 4 -^ 



«f>z- 



a 



(2) -2 



+ «3 



( 2), -3 



4 ” Cl 



( 2 ) -4 



I- b 






The difl'erence equation is 



y{n) = a^\x{n) + af\x{n — 1) + (t 2 ^lv(« — 2) + — 3) 

4- (74‘^.y(« — 4) 4- bf\\n — 3) 4- bf\v{n — 4). 



(3.43) 



(3.44) 



E. AUGMENTED HR NOTCH FILTER: P = 3 

From equations (3.12) and (3.13) in section A, we let p = 3. Those equations be- 
come: 



_i («o + ’ + « 2 - ")(1 + ^ 

Ji [z ) = 



(1 _ b,z-' _ z>2z-^)(l + aVT-‘ -F 



-2 



0),-' I J3) -2 



. 0 ) - 3 n 






1 - bf^z-^ - bf^z-^ 






T(z) 4^^ + 



T(z) 



1 _ bf>z-^ _ z,f ),-5 



(3.46) 



(3.47) 



(3.48) 

(3.49) 
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Manipulating (3.49) will result in the dilTerence equation, (3.53). 




(3.50) 



Y{z) — bf'’z~^Y{z) — bf^z ^y(z) = ao^^X(z) + af *A'(z) + af^z ^X{z) 

+ af ^z^^A^z) + ap^z~^A'(z) 



(3.51) 



+ A'(z) 



.< 3),-5 V, 



Y{z) = aQ\\\z) + af ^z ^X{z) + af^z ^X{z) + af^z ^X{z) + a^^z ^A'(z) 



y{n) = aP.x{n) + af \v(n — 1) + a 2 ^lv(« — 2) + af^x{/i — 3) + af^x{n — 4) 
+ af lv(« - 5) + bf^y{n - 4) + bf^y{n - 5). 



This realization is shown in Figure 6 on page 19. 

From Tables 1, 2 and 3 [Ref. 1; pp. 27S-279], we have 



ctf = ^f + Z,2 = -0.8712462 
af = b^^+ 2b^b2 = -0.61081465 
a^ = a, = 0.99 

^ = —0.0032993 

af^ = ao{bl + Z?2) + a^b^ + Oj = 0.01861246 
af^> = aoibl + Ib^bj) + a^{b^ + ^> 2 ) + «2^i = 0.00937445 
= a^{bi + Ib^bj) + a.ib^ + b^) = -0.66100766 
af> = a^_{b] + 2b^b^) = -0.6047065 
+ l>b\b2 + b\ = 0.65238232 
bf^ = b]b2 + Ib^bl = 0.59865944 . 

The transfer function is 



+ af z"^A'(z) + bf\-^ F(z) + bf\-^ Y{z) 



(3.52) 



c;f) = ^j =0.32993 




(3.55) 



The difference equation is 
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Figure 6. The Augmented HR Notch Filter: p = 3 
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(3.56) 



y{n) = OQ\x{n) + ~ I) + — 2) 

+ af lv(« — 3) + c[4^^j£-(« — 4) 

+ af\x(n - 5) + bf\{n - 4) 

+ bf\in - 5). 

F. STABILITY TESTING 

From section A if we start with a stable filter, we would like to be assured that the 
poles of the augmented-order filter all be outside the unit circle in the z~' plane. We de- 
sired this because, even though the added poles are cancelled by the zeroes of 
(1 4- 4- ••• 4- realization imperfections will prevent exact cancellation and the 

augmented filter would be unstable. In order to examine the stability question, we will 
use Jury's stability test as applied to successively higher-degree denominator polynomials 
in z as p increases [Ref 1: p. 281]. The results from this paper are summarized below: 

1. Condition for Stability of the Original HR Notch Filter: p = 0 

The original HR filter is described by equation (3.14) and to test for stability, 
we analyze the following characteristic polynomial resulting from it: 

F{z) = z- -b^z-b 2 = 0 (3.57) 

The conditions for stability are: 

1 — Z?] — ^2 > 0 ... ® 

1 +^i -Z)2>0... © (3.58) 

|/?2l < I ■•• © 

These conditions are shown graphically in Figure 7 on page 21. 

Circled numbers in the figure refer to the conditions of the same number in equation 
(3.58) 

Checking for stability of our example HR notch filter; we have from © of (3.58) 

1 — — /?2 > 0 
1 -(0.32993) -(-0.9801) >0 
1.65017 >0 



from © of (3.58) 
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Figure 7. Condition for Stability of the Original Filter [Ref. 1: p. 283] 



1 + 6 ] — > 0 
1 +(0.32993) -(-0.9801) >0 
2.31003 >0 



from © of (3.58) 

| 62 l<l 
10.98011 < 1 
0.9801 < 1 

Thus the conditions are all satisfied and the original filter is stable. 

2. Condition for Stability of the Augmented HR Notch Filter: p= 1 

The augmented HR notch filter, p= 1, is described by equation (3.33) and to test 
for stability, we analyze the following characteristic polynomial resulting from it: 



F{z) = 2 + 6, = 0 



(3.59) 



The conditions for stability are: 
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1 + > 0 ... © 
1 -Z>, >0... © 



(3.60) 



These conditions are shown graphically in Figure 8 on page 23. 

Circled numbers in the figure refer to the conditions of the same number in equation 
(3.60) 

Checking for stability of our example HR notch filter; we have from © of (3. 60) 

1 + > 0 
1 + 0.32993 >0 
1.32993 >0 

from © of (3.60) 

1 - > 0 
0.67007 > 0 

Thus the conditions are all satisfied and the augmented HR notch filter, p = 1, is stable. 
3. Condition for Stability of the Augmented HR Notch Filter: p = 2 

The augmented HR notch filter, p = 2, is described by equation (3.44) and to test 
for stability, we analyze the following characteristic polynomial resulting from it: 

F{z) = z'^ + b,z + {b] + b,) (3.61) 

The conditions for stability are: 

1 + + [ b ] + 162) > 0 ... © 

1 -Zj, +(/?,^ + ^2)>0- © (3-62) 

I + ^2 I < 1 ... ® 

These conditions are shown graphically in Figure 9 on page 24. 

Circled numbers in the figure refer to the conditions of the same number in equation 
(3.62) 

Checking for stability of our example HR notch filter; we have from © of (3. 62) 

1 + ^] + (^j + b-^ > 0 

1 + (0.32993) + [(0.32993)^ + ( -0.9801)] > 0 

0.4586838 >0 
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Figure 8. Condition for Stability of the Augmented Filter: p = 1 [Ref. 1: p. 284] 
from © of (3.62) 



1 - + [b] + > 0 

1 - (0.32993) + [(0.32993)^ + ( -0.9801)] < 0 

-0.201 1762 <0 



therefore the filter is unstable, from © of (3.62) 

\b\ + <1 

I (0.32993)^ + (-0.9801) I < 1 
0.8712462 < 1 

Note that © is violated and here the realization will be unstable. 

4. Condition for Stability of the Augmented HR Notch Filter: p = 3 

The augmented HR notch filter, p= 3, is described by equation (3.55) and to test 
for stability, we analyze the following characteristic polynomial resulting from it: 

F{z) = z V b,z^ + (bf + b 2 )z + (^^,^ + 2bib2) (3.63) 
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Figure 9. Condition for Stability of the Augmented Filter: p = 2 [Ref. 1: p. 285] 
The conditions for stability are; 



l+b^+ {b\ + b^) + {b\ + 26,^2) > 0 ... © 

\-b,+{b\-T b^) - {bi + 2b, bj) > 0 ... © 

\b,+2b,b2\ < 1 ... © 

1 (^f + 2b, bf - 1 1 > I + 2b,b2) - [b] + b^) I ... © 

These conditions are shown graphically in Figure 10 on page 25. 

Circled numbers in the figure refer to the conditions of the same number in equation 
(3.64) 

Checking for stability of our example HR notch filter; we have from ® of (3. 64) 

1 + /?] + (bf + ^> 2 ) "b (^1 "b ^ ^ 

1 + 0.32993 + [(0.32993)^ + ( -0.9801)] + [(0.32993)^ + 2(0.32993)( -0.9801)] > 0 

-0.15213085 <0 
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Figure 10. Condition for Stability of the Augmented Filter: p = 3 [Ref. 1: p. 2S6] 
therefore the filter is unstable, from © of (3.64) 



1 — fcj + [b] + ~ [b] + 2^1 ^’ 2 ) > 0 

1 - 0.32993 + [(0.32993)^ + ( -0.9S01)] - [(0.32993)^ + 2(0.32993)( -0.9801)] > 0 

0.40963S46 > 0 



from © of (3.64) 



1 + 2^,62! < 1 
I (0.32993)^ + 2(0.32993)( -0.9801)1 < 1 

0.61081465 < 1 



from © of (3.64) 



'\{b] + 2b^bj^ -\ \ > \b^{b]^2b,b:^)-{b\ + b:^)\ 

I {(0.32993)^ + 2(0.32993)( -0.9801)}^ - 1 1 < 

I (0.32993){(0.32993)^+ 2(0.32993)( -0.9801)} - {(0.32993)^ + ( -0.9801)} I 

0.62690546 <0.66972012 
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therefore the filter is unstable. 

Note that (D and ® are violated and therefore the realization will be unstable. 



G. FREQUENCY RESPONSE 

From the transfer functions (3.22), (3.33), (3.44), and (3.55), we use the program in 
Appendix A to calculate the magnitude and phase of the frequency response of the se- 
veral realizations. We plot the magnitude of these responses for the original filter and for 
the pipeline augmented that all of the frequency responses are the same from those fig- 
ures. 



1. Second Order HR Notch Filter: p = 0 

The magnitude response is shown in Figure 1 1 on page 27. 

2. The Augmented IIR Notch Filter: p = 1 

The magnitude response is shown in Figure 12 on page 28. 

3. The Augmented IIR Notch Filter: p = 2 

The magnitude response is shown in Figure 13 on page 29. 

4. The Augmented HR Notch Filter: p = 3 

The magnitude response is shown in Figure 14 on page 30. 

We know from section F that the example augmented HR notch filter for p = 
2, and 3 are unstable but we cannot see that from these magnitude responses. We have 
to find another way to show which one is stable or unstable. This will be shown in sec- 
tion H. 



H. UNIT SAMPLE RESPONSE 

From the difference equations (3.23), (3.34), (3.45), and (3.56), we use the program 
in Appendix B to calculate the output for a unit sample input to each of the pipeline 
realizations. We plot the outputs in Figures 15, 16, 17 and 18. We can see that the out- 
puts for Figure 15 (p = 0) and Figure 16 (p= 1) are stable, however the outputs for Fig- 
ure 17 (p= 2) and Figure 18 (p= 3) are unstable as predicted by the stability testing in 
section F. 

1. Second Order HR Notch Filter: p = 0 

The output is shown in Figure 15 on page 31. 

2. The Augmented HR Notch Filter: p = 1 

The output is shown in Figure 16 on page 32. 

3. The Augmented HR Notch Filter: p = 2 

The output is shown in Figure 17 on page 33. 
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Figure 11. Magnitude Response of Second Order HR Notch Filter: p = 0 
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MAGNITUDE RESPONSE 




Figure 12. Magnitude Response of the Augmented HR Notch Filter: p = 1 
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THETA (HAAS) 




Figure 13. Magnitude Response of the Augmented HR Notch Filter: p = 2 
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Figure 14. Magnitude Response of the Augmented HR Notch Filter: p = 3 
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OUTPUT SEQUENCE VCns) 




Figure 15. Output Sequence of Second Order HR Notch Filter: p = 0 
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SAHFLE « (ns) 




Figure 16. Output Sequence of The Augmented UR Notch Filter: p = 1 
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Figure 17. ‘ Output Sequence of The Augmented HR Notch Filter: p ~ 2 
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4. The Augmented HR Notch Filter: p = 3 

The output is shown in Figure 18 on page 35. 

I. PIPELINE IMPLEMENTATION 

All of the II R filters in sections B, C, D, and E can be developed into high-speed HR 
filters by applying pipeline techniques. We use Ka-stage pipeline add units and Km-stage 
pipeline multipliers to increase the speed of the filters. The Km stages are involved in the 
coefficient multiplies and the Ka stages are involved in the adds. The increase in order 
of the filter is P = 2Ka + Km — 2 [Ref 1: p. 280]. 

1. Second Order HR Notch Filter: p = 0 

To achieve a p = 0 realization, we cannot use pipeline modules with longer de- 
lays than Ka = 1 and Km = 0. 



P = 2Ka + Km -2 
= 2(l) + 0- 2 
= 0 

We have a realization of the original HR notch filter but with an additional 3 
delays. The equation is 

z~^Y=z~\aQX-\-a^z~'X+a2Z~^X+ b^z~^Y+ b2Z~^Y~\. (3.65) 

This realization is shown in Figure 19 on page 36. 

2. The Augmented HR Notch Filter: p = 1 

To achieve a p = 1 realization, we cannot use pipeline modules with longer de- 
lays than Ka = 1 and Km = 1 . 



P = 2Ka -I- Km - 2 
= 2 ( 1 )+ 1-2 
= 1 

We have a realization of the augmented HR notch filter, p = 1, but with an 
additional 4 delays. The equation is 

z-^Y=z~\aQX+a^z~^X+ a2Z~^X+ a2Z~^ X + b2Z~^ Y + £>32"^ T]. (3.66) 

This realization is shown in Figure 20 on page 37. 
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Figure 18. Output Sequence of The Augmented HR Notch Filter: p = 3 
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X 




Figure 19. Second Order HR Notch Filter: p = 0 uith Ka= 1, Km = 0 

3. The Augmented HR Notch Filter: p = 2 

To achieve a p=2 realization, we cannot use pipeline modules with longer de- 
lays than Ka = 1 and Km = 2. 



P = IKa + Km - 2 

= 2 ( 1 ) + 2-2 

_ T 


We have a realization of the augmented HR notch filter, p = 2, but with an 
additional 5 delays. The equation is 

2-5 Y = z“^[ooA' + a^2~^X + + a^z~^X + a 42 ~^A' + b^2~^ Y + b^2~^ T]. (3.67) 

This realization is shown in Figure 21 on page 38. 

4. The Augmented HR Notch Filter: p = 3 

To achieve a p=3 realization, we cannot use pipeline modules with longer de- 
lays than Ka = 1 and Km = 3. 
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Figure 20. The Augmented IIR Notch Filter: p = 1 with Ka= 1, Km= 1 
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Figure 21. Tlie Augmented IIR Notch Filter: p = 2 uith Ka= 1, Km=2 
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P = 2Ka + Km - 2 
= 2(l) + 3-2 
= 3 

We have a realization of the augmented HR notch filter, p = 3, but with an 
additional 6 delays. The equation is 

-ty -6 + a^z~^ X + 022~'^ X + X 

+ OiZ~^X + a^z~^X H- b^z~^ Y + b^z~ 

This realization is shown in Figure 22 on page 40. 



Table 4. LOOMIS/SINHA PIPELINE RECURSIVE DESIGN 



p-order 


Ka-add units 


Km-multipliers 


D (output delayed) 


Figure 


0 


1 


0 


3 


19 


1 


1 


1 


4 


20 


2 


1 




5 


21 




1 


J) 


6 


22 



(3.68) 



J. SUMMARY 

This chapter investigated a method for applying pipeline techniques to the design 
of a high-speed HR no' 'h filter. Using the Loomis Sinha method, HR notch filters op- 
erating at rates hitherto impossible can be designed. The general structure of the HR 
notch filter and the ^^’cthod of calculating the multiplier coefficients have been presented. 
The stability of the resulting realization has been investigated and a technique for as- 
certaining the stability of the realization has been presented. The delay and complexity 
of each realization are summarized in Table 4. 
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Figure 22. The Augmented HR Notch Filter: p = 3 >uth Ka= 1, Km= 3 
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IV. GNANASEKARAN'S STATE SPACE REALIZATIONS 



HR filter realizations are possible and most of them can be presented within the 
general time domain representation of the system under consideration. In continuous 
time it was convenient to represent a linear constant coefficient differential equation by 
a system of first-order differential equations for analysis purposes. It is, therefore, not 
surprising that it is useful to represent a linear constant coefficient difference equation 
by a system of first-order linear constant coefficient difference equations. Such a for- 
mulation results in a special form of representation of the difference equation called a 
state variable representation. 

A. STATE SPACE REPRESENTATION OF A NTH ORDER RECURSIVE 
DIGITAL FILTER 

We are used to working with a difference equation of the form 

;•(/;) = b^y{n - 1) -I- ••• -f b^v{n - X) + aQ.x{n) -k a^x{n - 1) -f ••• -f Oi^x{n - .1/) (4.1) 

or in the notation 



that is called the Nth-order form. An alternative is to describe an LTl system with N 
first-order difference equations and an output equation. This alternative form for de- 
scribing a linear time-invariant system is known as the state-space description and can 
be written as a matrix equation: 



In Eqs. (4.3) and (4.4), \{n) is a vector with A' rows and 1 column (an A'x 1 vector or 
column matri.x) of newly defined variables called the states, .v(«) is the system input, and 
y(«) its output. In the matrix state equation of Eq. (4.3) the matrix A must be an 




\{n) = Av(rt -1)4- B.v(/7) {state equation) 



(4.3) 



with the system output given by 



y(rt) = C^v(rt) 4- dx{n) {output equation). 



(4.4) 
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.V X ;V matrix and the B matrix an A’ x 1 vector. The scalar output >(«) is defined by Eq. 
(4.4) where C is a N x 1 column matrix and d is simply a real number that multiplies (or 
scales) the input .r(«). The elements in A, B, C, and d are constants that depend only 
upon the system parameters. 

B. WRITING THE STATE SPACE EQUATIONS 

To illustrate this approach to modeling of LTI systems we return to a version of the 
Nth-order difference equation described as a second order IIR filter, namely 

y{n) = b^y{n — 1) + b-^-{n — 2) + <3o^(/t) + — 1) + a 2 x{n — 2) (4.5) 

Y=b^2~^Y+b2Z~^Y+ QqX+ a^z~^X + a2Z~^X. (4.6) 

We define the state varible as: E,, V^. Then, we can establish the following re- 
lationships: 

}' = Y + b2Z~^ Y + aoX 4- a,z"’A' -I- a2Z~^X 

Fi = z"' Y (4.7) 

z-^r. 




Therefore the state equation is: 

(/ = b^2~^ r, + Z>2Z"’ V 2 + %X-ra^z~'X+a2Z~-X 
» 2 — ^ ’1. 

We express it in the matrix state equation: 

v(/i) = Av(/r — 1 ) -I- B.v(«) 

V = Az"’V -I- BA' 



>.■ 




1 




' v ; 




r . -1 , “2-1 




= 




2 




+ 




1 — 
I 




.1 




1 

1 




0 



That means: 



(4.8) 



(4.9) 

(4.10) 

(4.11) 



V = 



^^1 



IN 



(4.12) 
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The output equation will be: 



j(a 7) = C^v(n) + dx{n) 



(4.15) 




2 



(4.16) 



Y=c^\\+c^V^_ + dX. 



(4.17) 



C. MODIFIED STATE SPACE EQUATION 

From the state equation 



\’{n) = Av(rt — 1) + B.Y(rt), 



(4. IS) 



the computation of v(n) requires the immediate past state vector v(n-l). This recursive 
nature severely limits the speed that can be achieved by direct implementation of the 
state equation. The speed limits can be overcome by modifying this state equation so 
that v(n) depends not on the immediate past vector, v(/; — 1), but on r(« — L) for some 
L > 1. This can be accomplished as follows: 

For L = 1 the state equation is 



\\n) = Av (/7 — 1) + Bjr(rt) 



(4.19) 



V = A2"'V + BA. 



(4.20) 



Delaying equation (4.19) by 1, 2, and 3 yields 



\{n — 1) = Av(« — 2) + Bv(« — 1) 



(4.21) 



v(rt — 2) = Av(n — 3) + B.v(n — 2) 



(4.22) 



v(rt — 3) = Av(/7 — 4) + Bx{n — 3). 



(4.23) 



For L = 2 putting (4.21) into (4.22), the state equation, results is 



v(/?) = A[Av(w — 2) 4- B.v(n — 1)D + B.v(rt) 



(4.24) 



v(/;) = A^\{n — 2) + AB.r(« — I) + B.v(«) 



(4.25) 



V = A^z"^ V + ABz~' A' + BA^ 



(4.26) 



For L = 3 putting (4.22) into (4.25), the state equation, results is 



v(/j) = A^[Av(rt - 3) + Bx(rt - 2): + ABx(rt - 1) + Bx(«) (4.27) 



v(«) = A^v (/7 — 3) + A“Bx(rt — 2) + ABx(rt — I) + Bx(/?) 



(4.28) 



V = A^z"^' + A^Bz"^Y+ ABz”’.Y+ BA. 



(4.29) 



For Z- = 4 putting (4.23) into (4.28), the state equation, results is 

v(«) = A^[Av(« - 4) + B.v(/; - 3)] + A‘B.r(« - 2) + ABx(« - I) + Bx(n) (4.30) 

v(rt) = A^•(/? - 4) + A^B.v(/? - 3) + A^Bx(« - 2) + AB.v(/7 - I) + B.v(//) (4.3 1) 



One of the benefits of the state space realization is that if the original filter (L= I) 
is stable, all higher order realizations are also stable, and furthermore, A^ has improved 
error properties over A. This advantage of the state space method is discussed in the 
paper by Lti, Lee, and Messerschmitt [Ref 4]. 

D. SECOND ORDER HR FILTER STATE SPACE REALIZATIONS 

From the state space equation of the second order HR filter in section B and the 
modified state space equation in section C, we can develop techniques for high speed 
implementation of second order HR filters for several values of L. This section shows 
how to apply these concepts for L = 1,2, 3. and 4: 



V = A'^z'^V + A^Bz"^V + A‘Bz"lV + ABz"' A + BA. 



(4.32) 



Equation (4.32) can be written in the following form for any L> 1: 




(4.33) 
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1. L = I 



Beginning with equation (4.20), 

V = A2"'v + BX 



and substituting for V, A, and B from section B yields; 



V 2 



+ 



b { 


T 

N 


1 

1 


0 




— 1 
1 



Aq + d]Z ^ + CI2Z 
0 



X. 



(4.34) 



(4.35) 



Beginning with the output equation, equation (4.17), and referring to Figure 23, 
we note that the outpin of the multiplier d is dz~^X. This results is a revised output 
equation of 

Y = c, F, +C 2 V 2 + dz~^X. (4.36) 

Substituting the state equation into the output equation, results in 

Y={b,c, + C2)z-']\ 

+ b2C\Z V 2 

-f (7 qC|A 

+ ajCjz”^V 

+ (^2^1 

r = (),,c + c2)2-'F 

+ b2C-^z~'^Y 

+ c3[qCj A (4. 38) 

+ ajC,z X 
+ (<^2^1 "b ^ 

To demonstrate equality with the Y of the original HR filter as expressed in equation 
(3.23) for p = 0, let c, = 1, Cj = 0, and d=^ . Then the output equation is: 

Y = b^z~'^ Y + b2Z~^ Y + OqX + a^z~^X + a2Z~'^X. (4-39) 



By substituting: 
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(4.40) 






bn = b 

Clf \ = Cl 



12 



•0 — 1 



= a 
a') = a 



12 



•2 — ‘'M3 



we can write this equation as: 

Y = + a■^2Z~^X + 

This realization is shown in Figure 23 on page 47. 

2. L = 2 

Beginning with the state equation (4.26), 

V = A^z'-V + ABz"’A'+ BA' 
and substituting for V, A, and B from section B yields: 



'Vx 




' b , 


b { 


' b , 


b { 


^-2 


' M ' 


r 




_1 


0 _ 


_1 


0 _ 




. M . 



^1 ^’2 
1 0 



-1 , -2' 
^0 ^1^ + CI22 



+ 



-1 



ciq + a^z + a^z 



0 

- 2 "] 



-1 



0 










b'x d" bn 


1 

rsi 


—2 
7 ^ 


_F 2 . 




bi 


b2 





^■'1 

Vn 



+ 



+ 



2 

+ <2jZ?]Z + CZ2 ^i2 

_L -1 _L -2 

ao + ajz"' + 



z-‘A 



1 

1 




1 

] 





b\ d" ^2 


1 

CM 


-2 


MV 


bx 


/?2 


z 


M'^2. 



+ 



_ I _2 2 ' 

^0 + (^0^1 + +{a\b-^ + a-^z + fl 2 ^iz 

-1 , -2 , -3 

aoz + fljZ + G22 



A'. 



(4.41) 



(4.42) 



(4.43) 



(4.44) 



(4.45) 
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Figure 23. Second Order HR Filter: L = 1 
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Beginning with the output equation, equation (4.17), and referring to Figure 24, 
we note that the output of the multiplier d is dz~^X. This results is a revised output 
equation of 

F = c, Fj + C 2 F 2 + dz~^X. (4.46) 



Substituting the state equation into the output equation, results in 

Y=i{b] + b2)c, + b,c;]z-^V, 

+ "b ^ 2 ^ 2 ]^ ^2 

+ q^CyX 

+ + ai)ci + 

+ (^2^1 ^2^2 

Y^i{b\^b^)c,+b,c{\z-~Y 
+ [i>,Z52^i + ^ 2 *^ 2 ]^”^^^ 

+ OQCyX 

+ "b ^ 0 ^ 2 ]-^ 

+ + ^2^1 + aiC2)z~‘A' 

+ {Cl2b^C^ + ^2^2 + 



(4.47) 



(4.48) 



To demonstrate equality with the Y of the augmented HR filter as expressed in equation 
(3.34) for p = 1, let c, = 1, Cj = 0, and d =0. Then the output equation is: 

r= (fe^ + b2)z-^Y+b,b22-^Y+aoX+iaob, + a,)z~\Y 

—0 (4.‘4yj 

+ {ci\b^ + af2)z A + cijbiZ A. 



By substituting: 



^? + ^2 = ^ii 

^ 1^2 “ ^12 

Aq = <2n 
*^0^1 + <^i ~ <^12 
a^b^ + 02 = «i3 
^2^1 — <^\i 



(4.50) 



we can write this equation as: 
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(4.51) 



Y = Y + 6,2z"^ Y + «ii J + ai2z"'X 
+ a,3z“^A' + a^^z~^X. 



This is the same as the equation (3.34) for p= 1 found on page 15. This realization is 
shown in Figure 24 on page 50. 

3. L = 3 

Beginning with the state equation (4.29), 



V = + ABz~^X + BX 



(4.52) 



and substituting for V, A, and B from section B yields: 



1 

s 

1 


1 

1 


1 

K 

1 





bx + b^ ^>1^2 
b- 



b 

+ 

+ 



1 <^2 
bx + ^2 ^1^2 

bx ^2 

bx ^2 
1 0 





'bx b^' 


-3 


r^'.i 




\ 0 _ 


2 


— 1 

1 



, -2n 

^0 ^2^ 



-2 rr 

2 A 



, -1 , -2‘ 
Qq Q.\Z + 0-2^ 



-1 



Qq + a ^2 + ^22 

0 



0 

-21 



2“^A 






(4.53) 



V-, 





'bx{bl + 2b2) b2{bl + b,} 


-3 


>r 




bx 4" ^2 ^1 ^2 


2 


1"2. 



+ 



+ 



+ 



a,x{bx + ^ 2 ) + ^hi^\ + ^ 2 )^ ' + + ^ 2 )^ 



- 2 ' 



ciQbx + «i^]Z ' + a2bxZ 



-2 



r^v 



_ I 2 

ax^bx + UxbxZ + a^bxz 

1 -1 1 -2 

Qq + <3)Z + «22 



-1 I -2' 

Oq + ajZ + «22 



— ^ I'' 

z A 



X 



(4.54) 
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Q 



iB 



£] 



— (:^>> — 



Figure 24. The Augmented HR Filter: L = 2 








'b^(h} + 2b2) b2{b] + b2) 


^-3 


'K 




b\ + i>2 i>i^2 




Ji. 



+ 



«0 + («0^1 + «l)2 
+ (^1^? 



-1 



“ I " %^2 “ I " ^ 1^1 ^ 2 )^ 

£ Z ,^>2 + <^ 2 ^\)^~^ + (^ 2^1 + « 2 ^ 2 )-^~^ 
OqZ~^ + (( 3 o^i + «i)z~^ + {ci\b\ + a-^2~^ + a2b^z~^ 



X. 



(4.55) 



Beginning with the output equation, equation (4.17), and referring to Figure 25, 
we note that output of the multiplier d is dz~^X. This result is a revised output equation 
of 



Y = c, F, + C 2 K 2 + dz~^X. (4.56) 

Substituting the state equation into the output equation, results in 

y = [(i>f + 2^,^2 )c, + {b\ + ^>2)C2]z“^K, 

+ \_{b\b2 + b2)c^ + b^b^c^^z ^ V2 
-h 

+ [(«o^ + + aoC2]‘~'A’ (4.57) 

+ \_{oQb\ + OQb2 + + < 22 ) 0 , + (<2o^>i + a-^)c2~\z~^X 

+ [(<3]^^ + a^b2 + a2b\)C]^ + {a-^b^ + 02 )^ 23 ^“^ A' 

+ j|(i22^f "t" ^ 2 ^ 2 )^! "i" ^lb\^2 4" ^3^ 

}'= [(6^ + 26,/j2)<-'i + (b^\ + ^ 2 )^ 23 ^”'^' 

+ \_{b^ib2 + bl)c^ + b^b2C2^z~^Y 

+ C(«0^1 + ai)c, + OqC2~\z~'^X (4.5S) 

+ [(«o^f + «o^2 + ^\b\ + + («o^i + «i)c23^”^^ 

+ \_{a^b\ + a^b2 + a2^i)<^i + {<^\b\ + <^2)^23^~^'^' 

+ [(^2^? + ‘^ 2 ^ 2 )^^] + ‘^ 2 ^ 1^2 + d~\z~^X. 

To demonstrate equality with the Y of the augmented HR filter as e.xpressed in equation 
(3.45) for p = 2, let q = 1 , q = 0. and d = 0. Then the output equation is: 



51 



Y = (/^f + 2b^b2)z~^ Y + {bib2 + bl)z~^ Y + a^X + («o^i + «i)z"'A' 

C^o(^? ^2) ^2!]^ ^ ^2) ^2^1 3^ (‘^••59) 

+ « 2 (^? + b2)z~^X. 



By substituting: 



b] + 2b^b2 = b^ ] 

^1^2 + bl = b ^2 

(30 = «ii 
ao^i + = cJ]2 

2 

Oo('5’i + ^ 2 ) + ^ 1^1 + ^^2 = ^13 
+ ^ 2 ) + ^2^1 = «14 

« 2 (^? + b2) = a, 5 



(4.60) 



we can write this equation as; 



Y = b^z~^Y + b,;z~^Y + a,,X + anz~^ X 



'n^ 



Ml- 



+ ^132 ^X + <2|42 + ^2)52 A. 



(4.61) 



This is the same as the equation (3.45) for p=2 found on page 17. This realization is 
shown in Figure 25 on page 53. 

4. L = 4 

Beginning with equation (4.32), 



V = A^2~S^ + A^Bz~^X + A^Bz~-X+ AB2"’A'+ B.V 



(4.62) 



and substituting for V, A, and B from section B yields: 



52 
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~b,(b^ + 2 ^ 2 ) 


b2{b^ + bj) 


^ 2 ’ 


-a 

-7 


>r 




b^ + b2 


NJ 

1 


-1 o_ 


Z 


42 . 



+ 



+ 



+ 



+ 



'b,[b] + 2b^) b^{b\ + b^) 



+ ^>2 

2 

+ i>2 ^1^2 



b\h 



-1 . “2“1 
Oq + «lZ + 



-3 rr 

2 A 



-1 , - 2 ' 
^Zq 4- C?j2 + (^22 



-1 



^7q + a|2 + 0.2^ 

0 



0 

-2' 






-1 



X 



-1 , -2n 

(Zq + fliZ + Q2Z 



X 



(4.63) 



'bt+3b^b2 + bj b^b2{b] + 2b2) 


-4 


>r 


b,{bl + 2b2) b2{b] + b2) 


2 





+ 



+ 



+ 



+ 



cL^b^{b^ + 2Z?2) + Q.^b^{b^ + 2Z?2)^ + ^ 2^1 (^i ^^ 2 )^ 

%{b\ + /Z 2 ) + + ^ 2 )^”^ 



^oi^i ^’ 2 ) ^ 2 )-^ 

I 2 

«o^i + Q’^b-^z + a2^iz 






~ 1 “2 " 
i?o^i "(■ (z^b^z 4" 02^1^ 

CIq + <Z[2 ' + Q2Z 



z"‘x 



«0 -V a^z ’ + 
0 



A' 



z~^X 



(4.64) 





+ 2b\b2 + ^2 + 2 ^ 2 )' 


-4 






Z7,(^>f + 262 ) Z> 2 (^f + ^ 2 ) 


2 


42 . 



+ 



a^ + (oq^i + «i)2~' + (oq^? + ^^2 + ^ 1^1 + ^ 2 )^"^ 
+ (flo^i^ + 2ao^iZ?2 + + ^ 1^2 + <^2^\)^~^ 

+ + 2(3,^,Z?2 + 02 ^ 2 )-^”^ + (^2^f + 2a2b^b2)z~^ 

^7q2 + "I" ^ 1 )^ (^0^1 “I” ^0^2 "I" ^ 2 )^ 

+ {a^b^ + fli^2 + + i<^2^\ + ^ 2 ^ 2 )^”^ 



X. 



(4.65) 
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Beginning with the output equation, equation {4.17), and referring to Figure 26, 
we note that tlie output of the multiplier d is dz~^X. This results is a revised output 
equation of 



Y=c^]\ + C2V:^ + d2~^X. 
Substituting the state equation into the output equation, results in 



(4.66) 



(4.67) 



F= [(^f + ':>b]b2 + b])c^ + {b\ + 2b^b2)c:^']z~^\\ 

+ \_{b\b2 + 2b^b\)c^ + {b^bj + b\)c^z 
+ %c^X 

+ + ‘*0^2 T ^1^1 + ^ 2)^1 + (^ 0^1 + 

+ [(«o^? + 2^0 b\b2 + ^ 1^1 4" ^ 1^2 4” ^2^i)*'i 4" (^ 0^1 4" 4" (^\b\ 4" 02 )^ 2 !]“ 

4- 4- la^b\b2 4- a^b] + «2^2)^i + T ^i^-2 4- ciib\)c2^z~^X 

4- [( 02 ^? 4- 2a2b^b-^c^ 4- (« 2 ^i^ 4- a 2 b^C 2 4- d~\z~^X 

Y= [(6f + 3b^b2 + bl)c^ + (Z?f + 2b^b2)c2']z~^Y 
4- [(/)fZ?2 + 26,^-)ci + ib^b2 4- bl)c2']z~^ Y 

+ ^-(XiX 

-f [((ij/)| 4- a,)c, + a^C2'Jlz~'X 

4- [(uq^? + ^0^2 + ^ 1^1 4- ^' 2)^1 + (^ 0^1 + (4.6S) 

4- [(oq/^? 4- 2aQb]h2 4- a^b] 4- a]^2 + ^2^i)^'i 4- {%h\ + ^ 0^2 + ^ 1^1 + <^2)‘^2D"~^'^' 

4- [(a,6f 4- 2a^b^b2 + ^ 2 ^? 4- «2^2)‘^i 4- {a^b] 4- 01^2 4- a2b^)c2']z~^X 
+ [(«2^? 4- 2a2b]b2)c-^ + {a2b] 4- ct2i'2)^2 + d~\z~^X. 



To demonstrate equality with the Y of the augmented HR filter as expressed in equation 
(3.56) for p = 3, let c, — 1, Cj = 0, and (3'= 0. Then the output equation is: 



Y = (^f + Zb\b 2 + b\)z~^ Y + {btb 2 + 2b^b2)z~^ Y + aoY + {a^^b^ + a^)z~'X 
4- \_%{b] 4- Z?2) 4- a|i) 4- a2~\z 
4- [«o(^f + 2^,£'2) 4- a-^{b] + £^ 2 ) 4- a2b{]^z~^X 
+ 4" 2^j£>2) 4" ^ 2(^1 "f ^2)[]'^ "h ^ 2(^1 T '^b'^b-^z X. 



■ 2. - 5 , 



-1 



(4.69) 



Bv substituting: 



+ l)b\b2 + ^2 = i’ll 
b\b2 + Ib^bl = bi2 

Aq = fl]] 

Ofjb^ + «! = 0)2 

+ i>2) + ^li’i + ^2 = ^13 
aQ{b] + 2b^b2) + a^{b^ + 62 ) + a2b^ = 0,4 
a^{b^ + 2i?ii>2) + + i’2) = «i5 

«2(i’f + 2i>ii?2) = «16 

we can write this equation as: 

Y= i’i, 2 "^y+ i>, 2 z"^y+ aiiX+a^2Z~'X 

+ ni 3 Z~^r+ a, 4 z“^V (-1-71) 

+ A + ^i^z 



This is the same as the equation (3.56) for p = 3 on page 20. This realization is shown 
in Figure 26 on page 57. 

From equation (4.33) we can see that the feedback portion of the realization 
must have a delay of L stages in any loop. The general realization of that portion is 
shown in Figure 27 on page 58. 

E. PIPELINED IMPLEMENTATION 

All of the HR filters in section D can be developed into high-speed HR filters by 
applying pipeline techniques. We use Ka-stage pipeline add units and Km-stage pipeline 
multipliers to increase the speed of the filters. The Km stages are involved in the coeffi- 
cient multiplies and the Ka stages are involved in the adds. The net result is a (Km + 
Ka)-stagc pipeline multiply-add unit. The equations for the filter outputs for difierent 
L's are given below: 

I. L = I 

The equation is 

^-(8A'a+2A'm),r_ ^-(^Ka+lKm) 



-i- a]jAf+ a^2Z~'^X -I- 0]3z”^.V 



Y 



(4.72) 



This realizaticrris shown in Figure 28 on page 59. 
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Figure 27. Heart of the Pipeline Recursive Filter Realization 



2. L = 2 

The equation is 



(9A'fl+2A'm) y _ (9Aa-r2A^m) 



4 - + UjjZ 'a + ^132 A + CI14Z X 



(4.73) 



This realization is shown in Figure 29 on page 60. 
3. L = 3 



The equation is 



^ — (lO Afl+2A^) y __ ^~{10 Afl+2 AVt?) 



+ a^^X + a,2z“'A'+ 0,32“^ A' 
+ ^142 ^X + ^^152 X 



(4.74) 



This realization is shown in Figure 30 on page 61. 
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Figure 30. The Augmented HR Filter: L = 3 uith Ka, Km 



4. L = 4 



The equation is 



,-(llA:<J+2A:/?i)y _ ^-(WKa+lKm) 



+ a]]A^+fl]2Z *A^+ flj3Z 
+ a^^z~^X + «j5Z~^A^ + a^(,z~^X 



(4.75) 



This realization is shown in Figure 31 on page 63. 

From analysis of the feed back loop portion of the state space realization, as 
shown in Figure 27 and equation (4.33), we have seen that the delay around any loop 
must be L. A general realization for the feedback portion of the pipeline case is shown 
in Figure 32 on page 64. 

Therefore we can see that 



L = 2Ka + Km. (4.76) 

The following examples show the analysis of the cases for L = 1, 2, 3 and 4. 

5. Example for L = 1: 
a. Ka = 0, Km = I 



L = 2Ka + Km 
= 2(0) + 1 
= 1 

The only integer, non-negative values for Ka and Km which solve equation (4.76) are 
Ka = 0. Km= 1. This means we cannot use pipeline modules with longer delays than 
Ka = 0. Km — 1 . 

This realizes the original HR filter for L = 1, but with an additional 2 de- 
lays and the equation is 

r = Z~\b^ i2~‘ Y + Z?,2Z“^ Y + a, , Z + «,2Z"’ A' -f a^3Z~^X]. (4.77) 

This realization is shown in Figure 33 on page 65. 

6. E.xample for L = 2: 

The only integer, non-negative values for Ka and Km which solve equation 
(4.76) are Ka = 0, Km = 2 or Ka = I, Km = 0. 
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Figure 32. Heart of More General Pipeline Filter Realization 



a. Ka = 0, Km = 2 



L = 2Ka + Km 
= 2 ( 0 ) + 2 
= 2 

For these values, we use pipeline modules with delays Ka = 0, Km = 2 . 

This realizes the augmented HR filter for L = 2, but with an additional 4 
delays. The equation is 

z-^y = 2-^[6„z“^F+ ^2z"^r+ a,2z”'A'+ a^^z~^X+a^^z~^X']. (4.78) 

This realization is shown in Figure 34 on page 66. 
b. Ka = I, Km = 0 



L = 2Ka + Km 
= 2(l) + 0 
= 2 



For these values, we us" pipeline modules with delays Ka=\, Km = 0 . 
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Figure 33. Second Order HR Filter: L = 1 ^^ith Ka=0, Km= 1 
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This realizes the augmented IIR filter for L = 2, but tvith an additional 9 
delays. The equation is 

^ } = 2 'T+i>j22 } + fi [|22 + U ]22 + fl )42 (4.79) 

This realization is shown in Figure 35 on page 68. 

7. E.xample for L = 3: 

The only integer, non-negative values for Ka and Km which solve equation 
(4.76) are Ka = 0, Km = 3 or Afn = 1, Km = 1. 
a. Ka = 0, Km = 3 



L = IKa + Km 
= 2(0) + 3 



For these values, we use pipeline modules with delays Ka = 0. Km = 3 . 

This realizes the augmented IIR filter for L = 3, but with an additional 6 
delays. The equation is 



. J ,2“^ }• + 6,22”^ r + a, ,A' + a, 22"' A’ 
+ a, 32 "'^A'+ n, 42 ~^A'+ a, 52 ”'^A' 



(4.80) 



This realization is shown in Figure 36 on page 69. 
b. Ka = /, Km = I 



L = IKa 4- Km 
= 2 ( 1 ) + 1 
= j 

For these values, we use pipeline modules with delays Ka= 1, Km = 1 . 

This realizes the augmented HR filter for L = 3, but with an additional 12 
delays. The equation is 

bn^~^y+ «in^' + gj. 

-f fl) 32 ~^A' + £Ji 42 ~^A' H- U| 52 ~^A' 




This realization is shown in Figure 3 / on page 70. 
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Figure 36.' The Augmented HR Filter: L = 3 with Ka=0, Kni = 3 
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8. Example for L = 4: 

The only integer, non-negative values for Ka and Km which solve equation 
(4.76) are Ka = 0, Km = 4 and Ka= Km = 2 or /fu = 2, Km = 0. 
a. Ka = 0, Km = 4 



L = 2Ka -1- Km 
= 2(0) + 4 
= 4 



For these values, we use pipeline modules with delays Ka = 0, Km = 4 . 

This realizes the augmented HR filter for L = 4, but with an additional 8 
delays. The equation is 



-8 



Y = 



-8 

Z 



^]]Z } -f T+ cjjjA -I- fljjz 

+ «i 32~^A'4- a,4Z~^A' + a,5Z~‘^vf -I- a^^z~^X 



(4.82) 



This realization is shown in Figure 38 on page 72. 
b. Ka = /, Km = 2 



L = 2Ka + Km 
= 2(l) + 2 
= 4 

For these values, we use pipeline modules with delays Ka= 1, Km = 2 . 

This realizes the augmented HR filter for L = 4 but with an additional 15 
delays. The equation is 



-15 



Y = z 



-15 



^][Z } + ^ F ^ll'^ F ^12^ ^ 

+ a^jz~"X + fl]4z“^A' + U]5Z~^A' + a^(,z~^X 



(4.83) 



This realization is shown in Figure 39 on page 73. 
c. Ka = 2, Km = 0 



L = 2Ka + Km 
= 2 ( 2 ) + 0 
= 4 

For these values, we use pipeline modules with delays Ka = 2, Km = 0 . 

This realizes the augmented HR filter for L = 4, but with an additional 22 
delays. The equation is 
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73 



(4.84) 



-22 

2 



r= 



-22 



b^iZ~^Y+ b^2^~^Y+a^^X+a^2^'^X 
+ ^132 ^X H- ci^^z X + ci \^2 X + ^152 X 



This realization is shown in Figure 40 on page 75. 



Table 5. STATE SPACE REALIZATION SUMMARY 



L 


Ka-add units 


Km-multipliers 


D( output delayed) 


Figure 


1 


0 


1 


2 


31 


2 


0 


2 


4 


32 


2 


1 


0 


9 


33 




0 


3 


6 


34 




1 


1 


12 


35 


4 


0 


4 


S 


36 


4 


I 


2 


15 


37 


4 


2 


0 


22 


3S 



F. SUMMARY 

This chapter develops an alternate method for high speed implementation of state 
space digital filters. The proposed method implements A^ instead of A as seen in 
equation (4.33). The proposed method allows arbitrary A in the implementation, and it 
can take full advantage of the optimum state space filters along with the benefit of A^ 
realization. The delay and complexity of each realization are summarized in Table 5. 
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V. CONCLUSIONS 



A. SUMMARY 

This thesis describes pipeline realizations of the IIR notch filter. The transfer 
function and the coefficients that are appropiate for HR notch filter have been pre- 
sented. Realization of a high speed HR notch filter using a pipeline technique by Loomis 
and Sinha has been presented. An alternate technique in using state space realization 
by R. Gnanasekaran has also been presented. 

B. COMPARISON 

Table 4 in Chapter 111 summarizes the Loomis, Sinha pipeline recursive design and 
Table 5 in Chapter O/ suinmarizes Gnanasekaran's state space realization. It was shown 
in Chapter HI section I that we have only one choice of Ka and Km for each value of 
p and in Chapter IV section E that we have 1 choice of Ka and Km for L= 1, 2 choices 
of Ka and Km for L = 2 and for L = 3, and 3 choices of Ka and Km for L = 4. 

In the Loomis/Sinha pipeline recursive design we note that for p= 2, Ka= I, Km= 2, 
output delay = 5. the filter is unstable. However, in Gnanaskaran's state space realization 
for L = 4, Ka= 1. Km= 2, output delay= 15, the filter is stable. Thus we see that in cases 
where Loomis, Sinha realization arc unstable, a state space realization always exists and 
is stable. Hotve\er. it will require more hardware and have a greater overall delay. All 
of these results are summarized in Table 6. In that table, the column labelled D refers 
to the overall delay and the columns labelled irMlt and #Add refer to the number of 
multipliers and number of adders as shown in the appropriate figure. 

Note that for a given value of Ka and Km, the number of adders and multipliers 
required for the state space realization is larger than that required for the scalar. This 
disadvantage is offset by the guaranteed stability of the state space realization. 
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Table 6. SUMMARY OF REALIZATIONS 





State Space Realization 


Scalar Realization 


Ka 


Km 


L 


D 


Fig 


=Mlt 


#Add 


P 


D 


Fig 


#Mlt 


#Add 


0 


1 


1 


2 


31 


13 


11 












0 


2 


2 


4 


32 


15 


13 












1 


0 


2 


9 


33 


15 


13 


0 


3 


19 


5 


5 


0 


3 


-> 


6 


34 


17 


15 












1 


1 


'y 


12 


35 


17 


15 


1 


4 


20 


6 


6 


0 


4 


4 


8 


36 


19 


17 












1 


2 


4 


15 


37 


19 


17 


2 


5 


21 


7 


7 


7 


0 


4 


22 


38 


19 


17 












1 


3 












‘y 


6 


22 


8 
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APPENDIX A. PROGRAM AND DATA OF FREQUENCY RESPONSE 
A. PROGRAM FOR FREQUENCY RESPONSE CALCULATION 

DIGFREQ. FOR VERSION: 2/03/88 



PURPOSE: THIS PROGRAM COMPUTES THE FREOUENCY RESPONSE OF 

DISCRETE SYSTEMS. THE PROGRAM CONSISTS OF A MAIN 
PROGRAM THAT CONTROLS THE INPUT/OUTPUT AND THE 
SUBROUTINES dfresp AND coeff. SUBROUTINE dfresp 
COMPUTES THE FREOUENCY RESPONSE OF EACH SYSTEM. 
SUBROUTINE coeff ALLOWS THE USER THE OPTION OF 
GENERATING THE FILTER COEFFICIENTS OF THE SYSTEMS 
TO BE ANALYZED BY WRITING THE APPROPRIATE EOUATIONS. 

IF THE USER ELECTS TO GENERATE THE COEFFICIENTS BY 
USING THE SUBROUTINE coeff, THE EOUATIONS MUST BE 
WRITTEN INTO THE SUBROUTINE USING STANDARD FORTRAN 77 
STATEME.NTS. THE COEFFICIENTS MUST BE STORED IN THE 
ARRAYS bC) AND c() WHICH CORRESPOND RESPECTIVELY TO THE 
NUMERATOR AND DENOMINATOR TERMS OF THE SYSTEM EOUATION. 
THE USER CAN SELECT ONE OF TWO OPERATING MODES: BATCH 
OR TEST. IN BATCH MODE THE AMOUNT OF INTERFACE WITH 
THE USER IS MINIMIZED AND IT IS ASSUMED THAT THE INPUT 
DATA HAS BEEN STORED IN THE DEFAULT FILE ’ DIGFREO- IN ' . 

IN TEST MODE THE USER IS PROMPTED FOR THE NAME OF THE 
INPUT FILE OR HAS THE OPTION TO PERFORM A TRIAL RUN BY 
USING THE INPUT DATA STORED IN THE FILE ' DIGFREQ. TST’ . 

IT IS RECOMMENDED THAT FIRST-TIME USERS SELECT THE TEST 
MODE AND MAKE A TRIAL RUN WITH THE PRESTORED INPUT DATA. 
THE TEST NODE ECHOES THE INPUT DATA ONTO THE MONITOR TO 
ALLOW VERIFICATION OF ITS ACCURACY. THIS PROGRAM WILL 
COMPUTE THE FREQUENCY RESPONSE OF UP TO THREE SYSTEMS. 
FOR EACH SYSTEM, THE USER HAS THE OPTION OF HAVING THE 
OUTPUT EXPRESSED IN DECIBELS (db). THE OUTPUT OF THIS 
PROGRAM IS STORED IN TABULAR FORM IN THE FILE 
'DIGFREQ. OUT' AND IN A FORM SUITABLE FOR PLOTTING 
IN THE FILE 'DIGFREQ.DAT'. 












THIS PROGRAM ASSIMES THAT EACH DISCRETE SYSTEM IS MODELED BY THE 
EQUATION: H(z) = num/den WHERE: 

num = b(0)*z**T. -! b( l)*z**(L-l) + ... + b(L-l)*z + b(L) 

den = c(0)*z**N + c( l)*z**(N-l) + ... + c(N-l)*z + c(N) 

L = A NON-NEGATIVE INTEGER, THE DEGREE OF THE NUMERATOR 
POLYNOMIAL. 

N = A NON-NEGATIVE INTEGER, THE DEGREE OF THE DENOMINATOR 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



POLYNOMIAL. 

b(0)...b(L) = REAL COEFFICIENTS OF THE NUMERATOR TERMS. 
c(0)...c(N) = REAL COEFFICIENTS OF THE DENOMINATOR TERMS. 

THE INPUT PARAMETERS SHOULD BE STORED IN A FILE NAMED 
'DIGFREQ. IN'. ALL OF THE READ STATEMENTS USED BY THIS PROGRAM 
REQUIRE FORMATTED INPUT. PARTICULAR ATTENTION SHOULD BE PAID 
TO THE FORMATS, ESPECIALLY THE USE OF THE DECIMAL POINT TO 
DENOTE 'REAL' NUMBERS. THE INPUT PARAMETERS REQUIRED BY THE 
PROGRAM ARE LISTED BELOW. 



c 

c 


NAME 


TYPE 


RANGE (ARRAYS) 




RESTRICTIONS 


c 


numsys 


INTEGER 






1 


<= numsys <= 3 


c 


L 


INTEGER 








0 <= L <= 128 


c 


N 


INTEGER 








0 <= N <= 128 


c 


dsorce 


CHARACTER 








'F' OR 'S' 


c 


yscal 


CHARACTER 






1 


STD' OR 'LOG' 


c 


thetaO 


REAL 










c 


dlthta 


REAL 










c 


numpts 


INTEGER 






L 


<= numpts <= 101 


c 


b() 


REAL 


0, 


1, 2, ... , L 




0 <= L <= 128 


c 

r 


cQ 


REAL 


0, 


1, 2, ... , N 




0 <= N <= 128 


u 

c 


WHERE: 












c 


numsys = 


THE NUMBER OF 


DISTINCT SYSTEMS H(z) 


TO 


BE ANALYZED. 



C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 



THIS INTEGER VALUE MUST OCCUR AT THE TOP OF THE INPUT 
FILE. IT DELINEATES THE NUMBER OF SYSTEMS TO BE READ BY 
THE PROGRAM A.ND ANALYZED. FOR EACH SYSTEM (1, ... , numsys) 
THE PARAMETERS BELOW MUST APPEAR IN THE INPUT FILE. 

L = AN I.NTEGER VALUE SPECIFYING THE DEGREE OF THE NUMERATOR 
POLYNOMIAL. 

N = AN INTEGER VALUE SPECIFYING THE DEGREE OF THE DENOMINATOR 
POLYNOMIAL. 

dsorce = THE CHARACTER STRING 'F' OR 'S' DENOTING WHETHER THE 

SYSTEM COEFFICIENTS ARE TO BE READ FROM THE INPUT FILE (F) 
OR GENERATED (S) THROUGH USE OF THE SUBROUTINE coeff. 

yscal = A CHARACTER STRING SPECIFYING THE DESIRED MAGNITUDE OPTION: 
'STD' WILL PRODUCE STANDARD MAGNITUDE OUTPUT; 

'LOG' WILL PRODUCE MAGNITUDE EXPRESSED IN DECIBELS (db). 

thetaO = THE STARTING VALUE OF THETA (RADS) AS IN Z=EXP( J*THETA) . 

dlthta = THE INCREMENT OF THETA (RADIANS). 

numpts = THE NUMBER OF FREQUENCY POINTS FOR WHICH THE OUTPUT IS 
TO BE COMPUTED. 

b() = THE NUMERATOR COEFFICIENTS IN ORDER b(0), b(l), ..., b(L). 

IF dsorce = 'F' IS SELECTED THEN THE USER MUST SUPPLY THE 
L+1 NUMERATOR COEFFICIENTS IN THE FILE. IF dsorce = 'S' 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



THEN THE USER HAS ELECTED TO GENERATE THE NUMERATOR 
COEFFICIENTS BY WRITING THE APPROPRIATE FORTRAN STATEMENTS 
IN THE SPACE PROVIDED IN SUBROUTINE coeff. IF THIS METHOD 
OF DATA GENERATION IS ELECTED THE PROGRAM MUST BE RECOMPILED 
BEFORE EXECUTION. 

c() = THE DENOMINATOR COEFFICIENTS IN ORDER c(0), c(l), c(N). 

IF dsorce = 'F' IS SELECTED THEN THE USER MUST SUPPLY THE 
N+1 DENOMINATOR COEFFICIENTS IN THE FILE. IF dsorce = 'S' 
THEN THE USER HAS ELECTED TO GENERATE THE DENOMINATOR 
COEFFICIENTS BY WRITING THE APPROPRIATE FORTRAN STATEMENTS 
IN THE SPACE PROVIDED IN SUBROUTINE coeff. IF THIS METHOD 
OF DATA GENERATION IS ELECTED THE PROGRAM MUST BE RECOMPILED 
BEFORE EXECUTION. 

NOTE: THE INPUT FORMAT STATEMENTS OCCUR IN THE MAIN PROGRAM 

FOLLOWING THE CAPTION: ******** INPUT FORMAT ********. 

THE FORM OF THE INPUT DATA FILE IS; 

LINE // ENTRIES FORMAT 



C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c* 

c 

c 

c 

c 



1 

2 

3 

NOTE ] 
NOTE ^ 
NOTE : 

WHERE: 



numsys 

L, N, dsorce, yscal 
dlthta , thetaO ,numpts 
b(k) , k=0, 1,. . . ,L 
c(k), k=0,l,... ,N 



il 

i3,tll,i3,t21,al,t31,a3 
2fl0. 0,i3 
6fl0. 0 
6fl0. 0 



NOTES 1. 



2 . 



NN = 1 + (L/6 ROUNDED DOWN TO THE NEXT SMALLER INTEGER). 

ND = 1 + (N/6 ROUNDED DOWN TO THE NEXT SMALLER INTEGER). 

THE NEXT NN LINES ARE ONLY REQUIRED IF dsorce = 'F'. IF 

dsorce = ’ S’ THEN THE USER HAS ELECTED TO GENERATE THE 

L+1 NUMERATOR COEFFICIENTS IN THE SUBROUTINE coeff. 

THE USER MUST PROVIDE THE APPROPRIATE FORTRAN STATEMENTS 
IN SUBROUTINE coeff TO GENERATE THE VALUES FOR b(). 

THE NTIXT ND LINES ARE ONLY REQUIRED IF dsorce = 'F*. IF 

dsorce = 'S' THEN THE USER HAS ELECTED TO GENERATE THE 

NtI DENOMINATOR COEFFICIENTS IN THE SUBROUTINE coeff. 

THE USER MUST PROVIDE THE APPROPRIATE FORTRAN STATEMENTS 
IN SUBROUTINE coeff TO GENERATE THE VALUES FOR c(). 

FOR numsys > 1 THE FORMAT OF LINES 2. . . IS REPEATED. 

THE FORMAT f 10. 0 USED FOR INPUT DATA PERMITS THE DECIMAL 
POINT TO BE PLACED ANYWHERE IN THE FIELD OF 10 COLUMNS 
AND ALSO ALLOWS THE EXPONENTIAL FORMAT TO BE USED (EG. 
3146. 2 = 3, 1462E+03). 



THE OUTPUT DATA CREATED BY THE PROGRAM IS STORED IN TABULAR FORM 
IN THE FILE ' DIGFREQ. OUT' . ADDITIONALLY, THE OUTPUT DATA IS 
WRITTEN INTO THE FILE 'DIGFREQ.DAT' TO FACILITATE PLOTTING BY 



3. 

4. 
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ooooooooooooooonooooooooooooooooooooooooooooooo 



C A SEPARATE, USER SUPPLIED PROGRAM. THE FORMAT OF THE DATA IN 
C 'DIGFREQ.DAT' IS: el2.6, 2.x, el2. 6. THE FIRST ENTRY CORRESPONDS 

C TO THE ORDINATE VALUE (THETA) AND THE SECOND ENTRY THE ABSCISSA 
C VALUE (MAGNITUDE OR PHASE). ADDITIONAL HEADER INFORMATION IS 
C WRITTEN INTO THE DATA FILE TO ALLOW FOR CONTROL AND LABELING OF 
C EACH PLOT. 

C 

C 

Q******^v*******'***Vc*yc*******’*ycVf EXAMPLE **^%’y«’*yr***Vfyc**yc****ycVr***Vf*Vc’*Vc** 

THE INPUT PARAMETERS FOR THE SYSTEM DESCRIBED BELOW ARE STORED IN 
THE SAMPLE INPUT FILE ' DIGFREQ. TST' AND CAN BE USED FOR A TRIAL 
RUN IN THE TEST MODE. 



SYSTEM: H(z)=z/(z-0. 5) 

GOAL: TO OBTAIN THE FREQUENCY RESPONSE FOR THIS SYSTEM FROM 

THETA = 0.0 TO THETA = 3. 14159 (PI RADIANS) IN STEPS 
OF dlthta = PI/10. 0 

FOR THE SYSTEM DESCRIBED ABOVE THE INPUT FILE IS: 

1 

001 001 F STD 

0. 314159 0. 0 Oil 

1.0 0.0 

1.0 -0.5 



THE RESULTING OUTPUT DATA FILE: 'DIGFREQ. OUT' IS: 

INPUT DATA FOR SYSTEM # 1 

INPUT DATA SOURCEFILE: DIGFREQ. TST 
DEGREE OF NUMERATOR = 1 

DEGREE OF DENOMINATOR = 1 

dsorce = F 

NUMBER OF FREQUENCY POINTS = 11 MAGNITUDE OPTION = STD 

STARTING VALUE OF THETA = . OOOOOOE+00 

INCREMENT OF THETA = . 314159E+00 

THE NUMERATOR COEFFICIENTS b(0) ,b( 1). . . b(L) ARE 

. lOOOE+01 . OOOOE+00 



THE DENOMINATOR COEFFICIENTS c(0) ,c( 1). . . c(N) ARE 
. lOOOE+01 -.5000E+00 



OUTPUT DATA FOR SYSTEM // 1 



THETA MAGNITUDE PHASE 

(RADIANS) (DEGREES) 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



. OOOOOOE+00 
. 314159E+00 
. 628318E+00 
. 942477E+00 
. 125664E+01 
. 157080E+01 
. 188495E+01 
. 219911E+01 
. 251327E+01 
. 282743E+01 
. 314159E+01 



. 200000E+01 
. 182897E+01 
. 150588E+01 
. 122886E+01 
. 103088E+01 
. 894428E+00 
. 800894E+00 
. 737654E+00 
. 696900E+00 
. 674038E+00 
. 666667E+00 



. OOOOOOE+00 
164149E+02 
262677E+02 
29807E+02 
293546E+02 
265651E+02 
223862E+02 
173608E+02 
118186E+02 
597793E+01 
484184E-04 



END OF RUN, SYSTEM # 1 



FOR ILLUSTRATIVE PURPOSES THE COEFFICIENTS b( ) AND c( ) COULD 
HAVE BEEN GENERATED BY SPECIFYING dsorce = 'S' ANT) WRITING THE 
APPROPRIATE FORTRAN STATEMENTS INTO SUBROUTINE coeff. THE 
STATEMENTS THAT COULD BE USED TO ACCOMPLISH THIS ARE WRITTEN INTO 
THE SUBROUTINE BUT ARE 'COM^IENTED OUT'. 






;Wc^VV?^V^V^VVc Vc * V? V? Vf VfV? 



MAIN PROGRAM *^V*Vr‘)VVfV?^V*'>V')V';V'5V‘:V*‘;V‘>V*V(rVrVr‘)V*Vc';W-'3V 



character infile*12, mode*l, ylabl*13, dsorce*!, yscal*3 
real mh(lOl), ph(lOl), thetav(lOl), c(0:128), b(0: 128) 

C PROMPT USER FOR MODE: BATCH OR TEST. 

write(*, 1115) 
read(*, 1117) mode 

if ( (mode. eq. * ). or. (mode. eq. * y ^ ) ) then 

mode = ’ Y ’ 
write(*, 1118) 
read(*, 1119) infile 
else 

infiie = ’DIGFREQ. IN^ 
endif 

C UNIT=1 DEFINED AS INPUT FILE. UNITS=2,3 DEFINED AS OUTPUT FILES. 

open(unit=l , f ile=inf ile, status=* old* , iostat=ierr , err=999) 
open(unit=2,file=*DIGFREQ. OUT* ) 
open(unit=3,file=*DIGFREQ. DAT* ) 

C READ INPUT PARAMETERS AND CONDUCT ERROR CHECKS. 

read( 1 , 1000) numsys 
numplts = numsys* 2 
write( 3 , 2000) numplts 

if ( (numsys. It. 1) . or. (numsys. gt. 3) ) then 
write(*, 1122) numsys 

stop * Error, numsys must be in the range: 1 <= numsys <= 3. * 
endif 
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do 10 nsys=l, numsys 

data mh/101*o. 0/, ph/101*0. 0/, thetav/101*0. 0/ 
data b/129*0. 0/, c/129*0. 0/ 

read( 1,1001) L, N, dsorce, yscal 
read( 1,1002) dlthta, thetaO, numpts 

if((L. It. 0). or. (L. gt. 128)) then 
write(*, 1124) nsys , L 

stop * Error, L must be in the range: 0 <= L <= 128. ' 
elseif( (N. It. 0). or. (N. gt. 128)) then 
writeC*, 1125) nsys, N 

stop ’Error, N must be in the range: 0 <= N <= 128. * 
endif 

if ( (dsorce. eq. ’ F’ ). or. (dsorce. eq. ’ f ’ ) ) then 
dsorce = ’F’ 

elseif ( ( dsorce. eq. ’ S ’ ) . or. ( dsorce. eq. ’s’)) then 
dsorce = ’ S ’ 
else 

write(*, 1018) dcorce 

stop ’The allowed values for dsorce are: ’’s’’ or ’’F’’.’ 

endif 

if( (numpts. It. 1). or. (numpts. gt. 101)) then 
write(*, 1127) nsys, numpts 

stop ’Error, numpts must be in the range: 1 <= numpts <= 101. ’ 
endif 

if ( (yscal. eq. ’ STD’ ) . or. (yscal. eq. ’ std’ ) ) then 
yscal = ’STD' 
ylabl = ’ MAGNITUDE ’ 

elseif ( (vscal. eq. ’LOG’ ). or. (yscal. eq. ’ log’ )) then 
yscal = ’LOG’ 
ylabl = ’ MAGNITUDE (db)’ 
else 

write( - , 1128) yscal 

stop ’Error, yscal must be the string: ’’LOG’’ or ’’STD’’.’ 

endif 

C FOR dsorce = ’F’ READ THE COEFFICIENTS b( ) AND c( ) FROM THE INPUT 
C FILE. FOR dsorce = ’S’ CALL coeff TO GENERATE THE COEFFICIENTS. 

if (dsorce. eq. ’F’) then 

readd, 1003) (b(k),k=0,L) 
readd, 1003) (c(k),k=0,N) 
else 

call coef f (L,N,nsys ,b,c) 
endif 

C WRITE INPUT DATA INTO THE OUTPUT FILE: DIGFREQ. OUT. 

write( 2 , 1008) nsys 
write( 2 , 1010) infile 
write( 2 , 1110) L 
write( 2 , 1111) N 
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write(2, 1019) 
write( 2 , 1112) 
write(2, 1113) 
write( 2 , 1114) 
write( 2 , 1004) 
write( 2 , 1005) 
write( 2 , 1006) 
write( 2 , 1005) 
write( 2 , 1009) 
write( 2 , 1126) 
write( 2 , 1007) 



dsorce 

numpts, yscal 

thetaO 

dlthta 

(b(k),k=0,L) 

(c(k),k=0,N) 

nsys 

ylabl 



C FOR TEST MODE ECHO ALL INPUTS ONTO MONITOR (UNIT = *). 



if (mode. eq. * ) 

write(*, 1120) 

writer’s 1110) 

write(*,llll) 
write(*, 1019) 
write(*, 1112) 
writer’s 1113) 
write(*, 1114) 
write('% 1004) 
writeO’S 1005) 
write(* , 1006) 
write(* j 1005) 
write(*, 1123) 
pause ’END OF 
endif 



then 

nsys, infile 

L 

N 

dsorce 

numpts, yscal 

thetaO 

dlthta 

(b(k),k=0,L) 

(c(k),k=0,N) 

nsys 

RUN STRIKE <CR> WHEN READY TO CONTINUE. ’ 



C CALL dfresp TO CO''FUTE THE FREQUENCY RESPONSE. 

call df resp( b, c , mh,ph , L, N, thetaO , dlthta, thetav, numpts , yscal) 

C WRITE RESULTS INTO OUTPUT FILE: DIGFREQ.DAT. 

write( 3 , 2001) numpts 
write(3,*) ’MAGNITUDE RESPONSE’ 
write(3,*) ’THETA (RADS)’ 
write( 3, 2003) ylabl 
do 55 np=l, numpts 

write(3,2010) thetav(np), mh(np) 

55 continue 

write( 3 , 2001) numpts 
write(3,*) ’PHASE RESPONSE’ 
write(3,*) ’THETA (RADS)’ 
write(3,2003) ’ PHASE (DEG) ’ 
do 56 np=l, numpts 

write( 3 , 2010) thetav(np), ph(np) 

56 continue 



C WRITE RESULTS INTO OUTPUT FILE: DIGFREQ. OUT. 
do 150 np=l, numpts 
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write( 2 , 1013) thetav(np) , mh(np)j ph(np) 
150 continue 
write( 2, 1123) nsys 

10 continue 

write(*, 1121) 

999 close(unit=l) 
close(unit=2) 
close(unit=3) 

if ( ierr. gt. 0) then 
write(*, 1116) ierr 
endif 



Q**Vr**-jV***** INPUT FORMAT Vc^w»'‘.tvc*****^v****vc 



1000 formate il) 

1001 formate i3, t 11, i3,t21,al,t31,a3) 

1002 formate2fl0. 0,i3) 

1003 formate6fl0. 0) 



C 



V 1 V ^ V 



*V?’5Wf^Wc’‘>Wc'>Wc'5Wf^V**‘V^V****Vc* 



1004 format(t4,'THE NUMERATOR COEFFICIENTS b(0) ,b( 1). . . b(L) ARE: ',/) 

1005 format(6(2X,ell. 4) ,//) 

1006 0format(//,t4, 'THE DENOMINATOR COEFFICIENTS c( 0) , c( 1) . . . c(N) ' , 

1' ARE: ’ ,/) 

1007 format(t6, ' (RADIANS) ’ ,t38,' (DEGREES)' ,/) 

1008 format(tl6,' INPUT DATA FOR SYSTEM # ',il,//) 

1009 format(///,tl6, ' OUTPUT DATA FOR SYSTEM ’,il,/) 

1010 format(t4, 'INPUT DATA SOURCEFILE: ' ,al2) 

1013 format( t4, 3( el2. 6 ,4x) ) 

1018 format( lx, ' dsorce = ’ ,al,2x, 'Error, illegal value for dsorce. ' ) 

1019 format( t4 ,' dsorce = ',al) 

1110 format(t4, 'DEGREE OF NUMERATOR = ' ,i3) 

1111 format(t4, 'DEGREE OF DENOMINATOR = ',i3) 

1112 0format(t4, 'NUrSER OF FREQUENCY POINTS = ', i3 ,t39 ,' MAGNITUDE ' , 

1' OPTION = ' ,a3) 

1113 format(t4, 'STARTING VALUE OF THETA = ' ,el2. 6) 

1114 format(t4, 'INCREMENT OF THETA = ',el2. 6,/) 

1115 Oformatdx, 'DO YOU WISH TO RUN THIS PROGRAM IN TEST', 

1' MODE ? (Y/N) <CR> : ', ,) 

1116 Oformat(///, lx, 'ERROR OPENING INPUT FILE, PROGRAM TERMINATED.', 

1//, lx, 'ERROR CODE: ' ,i4, /////) 

1117 format(al) 

1118 Oformat(/////,lx, 'TYPE THE NAME OF YOUR DATA FILE FOLLOWED', 

1' BY <CR>. ' ,/,' IF YOU DESIRE TO MAKE A TEST RUN USING THE' , 

2' SAMPLE DATA ALREADY STORED',/,' IN THE FILE: DIGFREQ. TST' , 

3 ' TYPE: DIGFREQ. TST <CR>' ,/,' FILENAME: ', ,) 

1119 forraat(al2) 

1120 format(////,t4, 'SYSTEM #' ,il, ' INPUT DATA SOURCEFILE: ' ,al2) 

1121 0format(//,t4, 'TABULAR OUTPUT DATA IS STORED IN FILE: DIGFREQ. OUT' , 
l/,t4, 'PLOTTING DATA IS STORED IN FILE: DIGFREQ.DAT. ') 

1122 format( ////// ,t2 ,' The value of numsys is: ',il,'.') 
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1123 formate/, lx, 13( ' END OF RUN, SYSTEM# ' , il , 2x, 13( ' - ' ) , //) 

1124 0format(//////,t2, 'The degree(L) of the numerator for system ', 

1'# ’ ,il,' is : L = ' ,i3,'. ') 

1125 OformatC////// ,t2, 'The degree(N) of the denominator for system' 

1, ' ,# ' ,il, ' is : N = ' ,i3, ' . ' ) 

1126 formate ///,t8, 'THETA' , t21 ,al3 ,t40 ,' PHASE' ) 

1127 formate ////// ,t2 ,' The value of numpts for system ',il,' is: ' ,i3) 

1128 formate///// ,t2, 'The value of yscal is: ',a3,'.') 

2000 formate il) 

2001 formate i3) 

2003 formate al3) 

2010 formate el2. 6 , 2x, el2. 6) 

end 



C SUBROUTINE: dfresp 



C PURPOSE: 

C 

C 

C 

C 



THIS SUBROUTINE COMPUTES THE FREQUENCY RESPONSE OF 
THE SYSTEM. ALL FREQUENCY CALCULATIONS ARE IN RADIANS, 
HOWEVER THE OUTPUT IS CONVERTED TO DEGREES. 

THE OUTPUT FORMAT FOR EACH FREQUENCY INCREMENT IS: 
MAGNITUDEeM) PHASEeP) AS IN: M*EXPeJ*P). 



subroutine dfrespeb, c,mh,ph,L,N, theta0,dltht a, thetav, numpts , yscal) 

real mhenumpts), phenumpts), thetave numpts ) , imz, rez 

real be0:L), ce0:N) 

character yscal*3 

complex z, den, num, h, ci 

C DEFINE CONSTANTS. 

pi = 4 . 0 *atane 1 . 0) 
ci = ei. 0,0. 0) 

C ITERATE FROM thetaO, IN INCREMENTS OF dlthta. 

do 100 np=l , numpts 
num = ci*be0) 
den = ci*ce0) 

thetavenp) = thetaO + enp-l)*dlthta 
rez = cose thetavenp) ) 
imz = s ine thetave np) ) 
z = cmplxerez , imz) 

C CALCULATE NUMERA'^OR FOR GIVEN VALUE OF THETA, IF L > 0. 

if e L. gt. 0) then 
do 50 k=l, L 

num = z*num + ci'-bek) 

50 continue 

endif 
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o n o n 



C CALCULATE DENOMINATOR FOR GIVEN VALUE OF THETA, IF N > 0. 



if(N. gt. 0) then 
do 70 k=l, N 

den = z*den + ci*c(k) 

70 continue 

endif 

h = num/den 

C CONVERT COMPLEX VALUE 'h' INTO MAGNITUDE(mh) AND PHASE(ph) TERMS. 
C IF yscal = 'LOG' THEN CONVERT MAGNITUDE TO DECIBELS (db). 

C DIVIDE BY ZERO AVOIDED BY 'if' STATEMENTS. 

mh(np) = cabs(h) 
ifCyscal. eq. ' log' ) then 

if(mh(np). gt. 0. 00001) then 
mh(np) = 20. 0*logl0(mh(np) ) 
else 

mh(np) = -100.0 
endif 
endif 

if(abs(real(h)). It. 1. Oe-15) then 

if( abs( aimag(h) ). le. 1. Oe-15) ph(np)=0. 0 
if(aimag(h). gt. 1. Oe-15) ph(np)=90. 0 
if(aimag(h). It. -1. Oe-15) ph(np)=-90. 0 
else 

ph(np) = ( 180. 0/pi)*atan2(airaag(h) , real(h) ) 
endif 

100 continue 

return 

end 



C 



SUBROUTINE: coeff 



PURPOSE: THIS SUBROUTINE ALLOWS THE USER TO GENERATE THE 

NUMERATOR A.ND DENOMINATOR COEFFICIENTS THAT DESCRIBE 
EACH SYSTEM TO BE ANALYZED. IF dsorce = 'S' THEN 
THE MAIN PROGRAM WILL CALL THIS SUBROUTINE. 



subroutine coef f(L,N,nsys ,b,c) 
real b(0:L), c(0:N) 

pi = 4. 0*atan( 1. 0) 
el = L 
en = N 



C 

C 

C 






DEVELOP THE EQUATIONS TO GENERATE VALUES FOR THE ARRAYS b( ) AND c( ) 
IN THIS SPACE. THE STATEMENTS TYPED IN MUST FOLLOW STA.NDARD 
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FORTRAN 77 RULES AND MAY USE FORTRAN 77 INTRINSIC FUNCTIONS SUCH AS: 
SIN(), COSO, ABS()... AN EXAMPLE IS SHOWN BELOW. NOTE THAT THE 
VALUE nsys CAN BE USED TO DISTINGUISH BETWEEN SYSTEMS IF MORE THAN 
ONE SYSTEM (numsys > 1) IS TO BE ANALYZED. 

*** example *** 

if(nsys. eq. 1) then 
do 2 i=0, L 

b(i) = cos( i*pi/( 2. 0*el) ) 

2 continue 

do 3 i=0, N 

c( i) = cos(2. 0*i*pi/(3. 0*en)) 

3 cont inue 
endif 












return 

end 



B. MAGNITUDE, PHASE OF SECOND ORDER HR NOTCH FILTER: P = 0 

INPUT DATA FOR SYSTEM // 1 



INPUT DATA SOURCEFILE: digfreq. tst 
DEGREE OF NUMERATOR = 2 

DEGREE OF DENOMINATOR = 2 

dsorce = F 

NUMBER OF FREQUENCY POINTS = 100 MAGNITUDE OPTION = STD 
STARTING VALUE OF THETA = . OOOOOOE+00 

INCREMENT OF THETA = .280600E-01 

THE NUMERATOR COEFFICIENTS b( 0) ,b( 1). . . b(L) ARE: 

.9900E+00 -. 3299E+00 . 9900E+00 



THE DENOMINATOR COEFFICIENTS c(0) ,c( 1). . . c(N) ARE: 
. lOOOE+01 -.3299E+00 . 9801E+00 
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OUTPUT DATA FOR SYSTEM 1 



THETA 

(RADIANS) 

OOOOOOE+00 

280600E-01 

561200E-01 

841800E-01 

112240E+00 

140300E+00 

168360E+00 

196420E+00 

224480E+C0 

252540E+00 

280600E+00 

30S660E+00 

336720E+00 

364780E+00 

392840E+00 

420900E+00 

448960E+00 

477020E+00 

505080E+00 

533140E+00 

561200E+00 

589260E+00 

617320E+00 

645380E+00 

673440E+00 

701500E+00 

729560E+00 

757620E+00 

785680E+00 

S13740E+00 

841800E+00 

869860E+00 

897920E+00 

925980E+00 

954040E+00 

982100E+00 

101016E+01 

103822E+01 

106628E+01 

109434E+01 

112240E+01 

115046E+01 

117852E+01 

120658E+01 

123464E+01 

126270E+01 

129076E+01 



MAGNITUDE 



. 999939E+00 
. 999939E+00 
. 999939E+00 
. 999939E+00 
. 999938E+00 
. 999938E+00 
. 999937E+00 
. 999936E+00 
. 999935E+00 
. 999934E400 
. 999933E+00 
. 999931E+00 
. 999930E+00 
. 999928E+00 
. 999926E+00 
. 999923E+00 
. 999920E+00 
. 999917E+00 
. 999914E+00 
. 999910E+00 
. 999906E+00 
. 999902E+00 
. 999896E+00 
. 999891E+00 
. 999884E+00 
. 999876E+00 
. 999868E+00 
. 999858E+00 
. 999847E+00 
. 999835E+00 
. 999820E+00 
. 999S03E+00 
. 999783E+00 
. 999759E+00 
. 999731E+00 
. 999697E+00 
. 999655E+00 
. 999603E+00 
. 999537E+00 
. 999453E+00 
. 999343E+00 
. 999193E+00 
. 998985E+00 
. 998681E+00 
. 998213E+00 
. 997440E+00 
. 996024E+00 



PHASE 

(DEGREES) 

. OOOOOOE+00 
193947E-01 
388289E-01 
-.583436E-01 
779784E-01 
977754E-01 
117779E+00 
138030E+00 
158581E+00 
179474E+00 
200769E+00 
222516E+00 
244780E+00 
267620E+00 
291112E+00 
315332E+00 
340366E+00 
366300E+00 
393251E+00 
421327E+00 
450663E+00 
481405E+00 
513723E+00 
547804E+00 
583872E+00 
622182E+00 
663033E+00 
706764E+00 
753788E+00 
804587E+00 
859743E+00 
919964E+00 
986100E+00 
105921E+01 
114061E+01 
123199E+01 
133548E+01 
145389E+01 
159095E+01 
175176E+01 
194344E+01 
217626E+01 
246560E+01 
283559E+01 
332633E+01 
400969E+01 
502855E+01 
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131882E+01 

134688E+01 

137494E+01 

140300E+01 

143106E+01 

145912E+01 

148718E+01 

151524E+01 

154330E+01 

157136E+01 

159942E+01 

162748E+01 

165554E+01 

168360E+01 

171166E+01 

173972E+01 

176778E+01 

179584E+01 

182390E+01 

185196E+01 

188002E+01 

190808E+01 

193614E+01 

196420E+01 

199226E+01 

202032E+01 

204838E+01 

207644E+01 

210450E+01 

213256E+01 

216062E+01 

218868E+01 

221674E+01 

224480E+01 

227286E+01 

230092E+01 

232898E+01 

235704E+01 

238510E+01 

241316E+01 

244122E+01 

246928E+01 

249734E+01 

252540E+01 

255346E+01 

258152E+01 

260958E+01 

263764E+01 

266570E+01 

269376E+01 

272182E+01 

274988E+01 

277794E+01 



. 992994E+00 
. 984530E+00 
. 942857E+00 
. 381903E-01 
. 939926E+00 
. 984111E+00 
. 992867E+00 
. 995971E+00 
. 997413E+00 
. 998198E+00 
. 998672E+00 
. 998979E+00 
. 999190E+00 
. 999341E+00 
. 999452E+00 
. 999537E+00 
. 999603E+00 
. 999655E+00 
. 999698E+00 
. 999732E+00 
. 999761E+00 
. 999785E+00 
. 999805E+00 
. 999822E+00 
. 999837E+00 
. 999850E+00 
. 999861E+00 
. 999871E+00 
. 999880E+00 
. 999888E+00 
. 999894E+00 
. 999901E+00 
. 999906E+00 
. 999911E+00 
. 999915E+00 
. 999919E+00 
. 999923E+00 
. 999926E+00 
. 999929E+00 
. 999932E+00 
. 999934E+00 
. 999937E+00 
. 999939E+00 
. 999941E+00 
. 999942E+00 
. 999944E+00 
. 999945E+00 
. 999947E+00 
. 999948E+00 
. 999949E+00 
. 999950E+00 
. 999951E+00 
. 999952E+00 



671241E+01 
100259E+02 
194059E+02 
877665E+02 
. 200014E+02 
. 102595E+02 
. 687161E+01 
. 516115E+01 
. 413001E+01 
. 344006E+01 
. 294546E+01 
. 257310E+01 
. 228226E+01 
. 204852E+01 
. 185629E+01 
. 169521E+01 
. 155808E+01 
. 143978E+01 
. 133653E+01 
. 124550E+01 
. 116455E+01 
. 109198E+01 
. 102647E+01 
. 966960E+00 
. 912579E+00 
. 862632E+00 
. 816531E+00 
. 773792E+00 
. 734015E+00 
. 696843E+00 
. 661987E+00 
. 629193E+00 
. 598246E+00 
. 568950E+00 
. 541144E+00 
. 514683E+00 
.489434E+00 
. 465285E+00 
. 442146E+00 
. 419918E+00 
. 398521E+00 
. 377889E+00 
. 357941E+00 
. 338639E+00 
. 319915E+00 
. 301724E+00 
. 284024E+00 
. 266771E+00 
. 249929E+00 
. 233459E+00 
. 217332E+00 
. 201519E+00 
. 185987E+00 



END OF RUN, SYSTEM # 1 
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C. MAGNITUDE, PHASE OF THE AUGMENTED HR NOTCH FILTER: P= 1 



INPUT DATA FOR SYSTEM # 1 



INPUT DATA SOURCEFILE: digfreq. tst 
DEGREE OF NUMERATOR = 3 

DEGREE OF DENOMINATOR = 3 

dsorce = F 

NUMBER OF FREQUENCY POINTS = 100 MAGNITUDE OPTION = STD 
STARTING VALUE OF THETA = . OOOOOOE+00 

INCREMENT OF THETA = . 280600E-01 

THE NUMERATOR CO: '^FICIENTS b( 0) ,b( 1) . . . b( L) ARE: 

.9900E+00 3299E-02 . 8811E+00 . 3266E+00 



THE DENOMINATOR COEFFICIENTS c(0) ,c( 1). . . c(N) ARE: 
. lOOOE+01 .OOOOE+00 . 8712E+00 . 3234E+00 



OUTPUT DATA FOR SYSTEM // 1 



THETA 

(RADIANS) 

. OOOOOOE+00 
. 280600E-01 
. 561200E-01 
. 841800E-01 
. 112240E+00 
. 140300E+00 
. 168360E+00 
. 196420E+00 
. 224480E+00 
. 252540E+00 
. 280600E+00 
. 308660E+00 
. 336720E+00 
. 364780E+00 
. 392840E+00 
. 420900E+00 
. 448960E+00 
. 477020E+00 



MAGNITUDE 



. 999940E+00 
. 999940E+00 
. 999939E+00 
. 999939E+00 
. 999939E+00 
. 999938E+00 
. 999937E+00 
. 999936E+00 
. 999935E+00 
. 999934E+00 
. 999933E+00 
. 999932E+00 
. 999930E+00 
. 999928E+00 
. 999926E+00 
. 999923E+00 
. 999921E+00 
. 999918E+00 



PHASE 

(DEGREES) 

. OOOOOOE+00 
193950E-01 
388293E-01 
583432E-01 
779798E-01 
977766E-01 
117780E+00 
138031E+00 
158579E+00 
179475E+00 
200768E+00 
222515E+00 
244780E+00 
267618E+00 
291112E+00 
315332E+00 
340364E+00 
366303E+00 



91 



505080E+00 

533140E+00 

561200E+00 

589260E+00 

617320E+00 

645380E+00 

673440E+00 

701500E+00 

729560E+00 

757620E+00 

785680E+00 

813740E+00 

841800E+00 

869860E+00 

897920E+00 

925980E+00 

954040E+00 

982100E+00 

1C1016E+01 

103822E+01 

106628E+01 

109434E+01 

112240E+01 

115046E+01 

117852E+01 

120658E+01 

123464E+01 

126270E+01 

129076E+01 

131SS2E+01 

134&88E+01 

137494E+01 

140300E+01 

143106E+01 

145912E+01 

148718E+01 

151524E+01 

154330E+01 

157136E+01 

159942E+01 

162748E+01 

165554E+01 

168360E+01 

171166E+01 

173972E+01 

176778E+01 

179584E+01 

182390E+01 

185196E+01 

188002E+01 

190808E+01 

193614E+01 

196420E+01 

199226E+01 

202032E+01 

204838E+01 



. 999914E+00 
. 999911E+00 
. 999906E+00 
. 999902E+00 
. 999897E+00 
. 999891E+00 
. 999884E+00 
. 999877E+00 
. 999868E+00 
. 999859E+00 
. 999848E+00 
. 999835E+00 
. 999820E+00 
. 999803E+00 
. 999783E+00 
. 999759E+00 
. 999731E+00 
. 999697E+00 
. 999655E+00 
. 999603E+00 
. 999538E+00 
. 999454E+00 
. 999343E+00 
. 999194E+00 
. 998985E+00 
. 99868 lE+00 
. 998214E+00 
. 997441E+00 
. 996025E+00 
. 992995E+00 
. 984532E+00 
. 942860E+00 
. 382032E-01 
. 939922E+00 
. 984110E+00 
. 992866E+00 
. 995970E+00 
. 997413E+00 
. 998197E+00 
. 99867 lE+00 
. 998979E+00 
. 999189E+00 
. 999340E+00 
. 999452E+00 
. 999537E+00 
. 999603E+00 
. 999655E+00 
. 999697E+00 
. 999732E+00 
. 999761E+00 
. 999785E+00 
. 999805E+00 
. 999822E+00 
. 999837E+00 
. 999850E+00 
. 999861E+00 



393250E+00 
421327E+00 
450665E+00 
481406E+00 
513721E+00 
-.547805E+00 
583875E+00 
622183E+00 
663032E+00 
706764E+00 
753788E+00 
804590E+00 
859746E+00 
919964E+00 
986101E+00 
105921E+01 
114061E+01 
123199E+01 
133548E+01 
145389E+01 
159095E+01 
175176E+01 
194344E+01 
217625E+01 
246560E+01 
283559E+01 
332632E+01 
400968E+01 
502853E+01 
671239E+01 
100259E+02 
194058E+02 
877594E+02 
. 200013E+02 
. 102595E+02 
. 687158E+01 
. 516112E+01 
. 412999E+01 
. 344004E+01 
. 294544E+01 
. 257309E+01 
. 228225E+01 
. 204850E+01 
. 185627E+01 
. 169520E+01 
. 155807E+01 
. 143977E+01 
. 133652E+01 
. 124549E+01 
. 116454E+01 
. 109197E+01 
. 102646E+01 
. 966947E+00 
. 912566E+00 
. 862621E+00 
. 816522E+00 
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. 207644E+01 
. 210450E+01 
. 213256E+01 
. 216062E+01 
. 218868E+01 
. 221674E+01 
. 224480E+01 
. 227286E+01 
. 230092E+01 
. 232898E+01 
. 235704E+01 
. 238510E+01 
. 241316E+01 
. 244122E+01 
. 246928E+01 
. 249734E+01 
. 252540E+01 
. 255346E+01 
. 258152E+01 
. 260958E+01 
. 263764E+01 
. 266570E+01 
. 269376E+01 
. 272182E+01 
. 274988E+01 
. 277794E+01 



. 999871E+00 
. 999880E+00 
. 999887E+00 
. 999894E+00 
. 999900E+00 
. 999906E+00 
. 999911E+00 
. 999915E+00 
. 999919E+00 
. 999923E+00 
. 999926E+00 
. 999929E+00 
. 999932E+00 
. 999934E+00 
. 999936E+00 
. 999939E+00 
. 999940E+00 
. 999942E+00 
. 999944E+00 
. 999945E+00 
. 999947E+00 
. 999948E+00 
. 999949E+00 
. 999950E+00 
. 999951E+00 
. 999952E+00 



. 773783E+00 
. 734003E+00 
. 696833E+00 
. 661978E+00 
. 629187E+00 
. 598235E+00 
. 568938E+00 
. 541135E+00 
. 514669E+00 
. 489423E+00 
. 465277E+00 
. 442138E+00 
. 419907E+00 
. 398514E+00 
. 377878E+00 
. 357937E+00 
. 338629E+00 
. 319905E+00 
. 301716E+00 
. 284016E+00 
. 266763E+00 
. 249920E+00 
. 233452E+00 
. 217326E+00 
. 201509E+00 
. 185980E+00 



END OF RUN, SYSTEM ^ 1 



D. MAGNITUDE, PHASE OF THE AUGMENTED HR NOTCH FILTER: P 

INPUT DATA FOR SYSTEM 1 



INPUT DATA SOURCr.FILE: digfreq. tst 
DEGREE OF NUMERAL OR = 4 

DEGREE OF DENOMINATOR = 4 

dsorce = F 

NUMBER OF FREOUENGY POINTS = 100 MAGNITUDE OPTION = STD 
STARTING VALUE OF THETA = . OOOOOOE+00 

INCREMEOT OF THETA = . 280600E-01 

THE NU’MEPvATOR COEFFICIENTS b(0) ,b( 1). . . b(L) ARE: 

.9900E+00 -.3299E-02 . 1861E-01 . 6141E+00 8625E+00 



THE DENOMINATOR COEFFICIENTS c( 0) , c( 1) . . . c(N) ARE: 

. lOOOE+01 . OOOOE+00 . OOOOE+00 . 6108E+00 8539E+00 
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OUTPUT DATA FOR SYSTEM # 1 



THETA 

(RADIANS) 

OOOOOOE+00 

280600E-01 

561200E-01 

841800E-01 

112240E+00 

140300E+00 

168360E+00 

196420E+00 

224480E+00 

252540E+00 

280600E+00 

308660E+00 

336720E+00 

364780E+00 

392840E+00 

420900E+00 

448960E+00 

477020E+00 

505080E+00 

533140E+00 

561200E+00 

589260E+00 

617320E+00 

645380E+00 

673440E+00 

701500E+00 

729560E+00 

757620E+00 

785680E+00 

813740E+00 

S41800E+00 

869860E+00 

897920E+00 

925980E+00 

954040E+00 

982100E+00 

101016E+01 

103822E+01 

106628E+01 

109434E+01 

112240E+01 

115046E+01 

117852E+01 

120658E+01 

123464E+01 



MAGNITUDE 



. 999940E+00 
. 999940E+00 
. 999940E+00 
. 999940E+00 
. 999939E+00 
. 999938E+00 
. 999938E+00 
. 999937E+00 
. 999936E+00 
. 999934E+00 
. 999933E+00 
. 999932E+00 
. 999930E+00 
. 999928E+00 
. 999926E+00 
. 999923E+00 
. 999920E+00 
. 999917E+00 
. 999914E+00 
. 999910E+00 
. 999906E+00 
. 999901E+00 
. 999896E+00 
. 999891E+00 
. 999884E+00 
. 999876E+00 
. 999868E+00 
. 999858E+00 
. 999847E+00 
. 999835E+00 
. 999820E+00 
. 999803E+00 
. 999783E+00 
. 999759E+00 
. 999731E+00 
. 999697E+00 
. 999655E+00 
. 999603E+00 
. 999537E+00 
. 999453E+00 
. 999343E+00 
. 999193E+00 
. 998985E+00 
. 998681E+00 
. 998213E+00 



PHASE 

(DEGREES) 

. OOOOOOE+00 
194012E-01 
388423E-01 
583613E-01 
780022E-01 
978035E-01 
117802E+00 
138056E+00 
158606E+00 
179500E+00 
200794E+00 
222545E+00 
244805E+00 
267638E+00 
291136E+00 
315351E+00 
340386E+00 
366321E+00 
393270E+00 
421345E+00 
450679E+00 
-.481422E+00 
513734E+00 
547815E+00 
583885E+00 
622197E+00 
663044E+00 
706773E+00 
753796E+00 
804595E+00 
859751E+00 
919967E+00 
986102E+00 
105921E+01 
114062E+01 
123199E+01 
133548E+01 
145388E+01 
159094E+01 
175175E+01 
194343E+01 
217625E+01 
246558E+01 
283557E+01 
332630E+01 
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126270E+01 

129076E+01 

131882E+01 

134688E+01 

137494E+01 

140300E+01 

143106E+01 

145912E+01 

148718E+01 

151524E+01 

154330E+01 

157136E+01 

159942E+01 

162748E+01 

165554E+01 

168360E+01 

171166E+01 

173972E+01 

176778E+01 

179584E+01 

182390E+01 

185196E+01 

188002E+01 

19080SE+01 

193614E+01 

196420E+01 

199226E+01 

202032E+01 

204838E+01 

207644E+01 

210450E+01 

213256E+01 

216062E+01 

218S68E+01 

221674E+01 

224480E+01 

227286E+01 

230092E+01 

232898E+01 

235704E+01 

238510E+01 

241316E+01 

244122E+01 

246928E+01 

249734E+01 

252540E+01 

255346E+01 

258152E+01 

260958E+01 

263764E+01 

266570E+01 

269376E+01 

272182E+01 

274988E+01 

277794E+01 



. 997440E+00 
. 996024E+00 
. 992994E+00 
. 984530E+00 
. 942856E+00 
. 381887E-01 
. 939924E+00 
. 984111E+00 
. 992867E+00 
. 995970E+00 
. 997413E+00 
. 998198E+00 
. 998671E+00 
. 998979E+00 
. 999189E+00 
. 999340E+00 
. 999452E+00 
. 999537E+00 
. 999603E+00 
. 999655E+00 
. 999697E+00 
. 999732E+00 
. 999761E+00 
. 999785E+00 
. 999805E+00 
. 999822E+00 
. 999837E+00 
. 999850E+00 
. 999861E+00 
. 999871E+00 
. 999880E+00 
. 999887E+00 
. 999894E+00 
. 999900E+00 
. 999906E+00 
. 999911E+00 
. 999915E+00 
. 999919E+00 
. 999923E+00 
. 999926E+00 
. 999929E+00 
. 999932E+00 
. 999934E+00 
. 999936E+00 
. 999938E+00 
. 999940E+00 
. 999942E+00 
. 999944E+00 
. 999945E+00 
. 999947E+00 
. 999948E+00 
. 999949E+00 
. 999950E+00 
. 999951E+00 
. 999952E+00 



400966E+01 
502849E+01 
671234E+01 
100258E+02 
194057E+02 
877335E+02 
. 200011E+02 
. 102594E+02 
. 687151E+01 
. 516108E+01 
.412995E+01 
. 34400 lE+01 
. 294542E+01 
. 257306E+01 
. 228223E+01 
. 204849E+01 
. 185626E+01 
. 169518E+01 
. 155806E+01 
. 143976E+01 
. 133650E+01 
. 124548E+01 
. 116454E+01 
. 109196E+01 
. 102646E+01 
. 966945E+00 
. 912560E+00 
. 862616E+00 
. 816516E+00 
. 773779E+00 
. 734001E+00 
. 696832E+00 
. 661977E+00 
. 629182E+00 
. 598231E+00 
. 568938E+00 
. 541134E+00 
. 514669E+00 
, 489424E+00 
. 465279E+00 
. 442136E+00 
, 419910E+00 
. 398511E+00 
. 377878E+00 
. 357935E+00 
. 338632E+00 
. 319908E+00 
. 301717E+00 
. 284017E+00 
. 266763E+00 
. 249920E+00 
. 233450E+00 
. 217325E+00 
. 201511E+00 
. 185983E+00 
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END OF RUN, SYSTEM // 1 



E. MAGNITUDE, PHASE OF THE AUGMENTED HR NOTCH FILTER: P = 3 



INPUT DATA SOURCEFILE: digfreq. tst 
DEGREE OF NUMERATOR = 5 

DEGREE OF DENOMINATOR = 5 

dsorce = F 

NUMBER OF FREQUENCY POINTS = 100 MAGNITUDE OPTION = STD 
STARTING VALUE OF THETA = . OOOOOOE+00 

INCREMENT OF THETA = . 280600E-01 

THE NUMERATOR COEFFICIENTS b( 0) ,b( 1). . . b(L) ARE; 

. 9900E+00 -.3299E-02 . 1861E-01 . 9374E-02 6610E+00 

6047E+00 



THE DENOMINATOR COEFFICIENTS c(0) ,c( 1). . . c(N) ARE: 

. lOOOE+01 . OOOOE+00 . OOOOE+00 . OOOOE+00 6524E+00 
5987E+00 



INPUT DATA FOR SYSTEM # 1 



OUTPUT DATA FOR SYSTEM // 1 



THETA 

(RADIANS) 



MAGNITUDE 



PHASE 

(DEGREES) 



. OOOOOOE+00 
. 280600E-01 
. 561200E-01 
. 841800E-01 
. 112240E+00 
. 140300E+00 
. 168360E+00 
. 196420E+00 
. 224480E+00 
. 252540E+00 



999937E+00 . OOOOOOE+00 
999938E+00 194494E-01 
999938E+00 388970E-01 
999938E+00 -.584018E-01 
999938E+00 780285E-01 
999938E+00 978157E-01 
999937E+00 117817E+00 
999936E+00 138065E+00 
999935E-f00 158610E+00 
999934E+00 179497E+00 
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280600E+00 

308660E+00 

335720E+00 

364780E+00 

392840E+00 

420900E+00 

448960E+00 

477020E+00 

505080E+00 

533140E+00 

561200E+00 

589260E+00 

617320E+00 

645380E+00 

673440E+00 

701500E+00 

729560E+00 

757620E+00 

785680E+00 

813740E+00 

841800E+00 

869860E+00 

897920E+00 

925980E+00 

954040E+00 

9S2100E+00 

101016E+01 

103822E+01 

106628E+01 

109434E+01 

112240E+01 

115046E+01 

117852E+01 

120658E+01 

123464E+01 

126270E+01 

129076E-i-01 

131882E+01 

134688E+01 

137494E+01 

140300E+01 

143106E+01 

145912E+01 

148718E+01 

151524E+01 

154330E+01 

157136E+01 

159942E+01 

162748E+01 

165554E+01 

168360E+01 

171166E+01 

173972E+01 

176778E+01 

179584E+01 

182390E+01 



. 999933E+00 
. 999931E+00 
. 999929E+00 
. 999927E+00 
. 999925E+00 
. 999923E+00 
. 999920E+00 
. 999917E+00 
. 999914E+00 
. 999910E+00 
. 999906E+00 
. 999902E+00 
. 999896E+00 
. 999891E+00 
. 999884E+00 
. 999876E+00 
. 999868E+00 
. 999858E+00 
. 999847E+00 
. 999835E+00 
. 999820E+00 
. 999803E+00 
. 999783E+00 
. 999759E+00 
. 999731E+00 
. 999697E+00 
. 999655E+00 
. 999603E+00 
. 999537E+00 
. 999453E+00 
. 999343E+00 
. 999193E+00 
. 998985E+00 
. 998681E+00 
. 998213E+00 
. 997440E+00 
. 996024E+00 
. 992994E+00 
. 984530E+00 
. 942856E+00 
. 381933E-01 
. 939924E+00 
. 984110E+00 
. 992867E+00 
. 995970E+00 
. 997413E+00 
. 998198E+00 
. 998671E+00 
. 998979E+00 
. 999189E+00 
. 999340E+00 
. 999452E+00 
. 999537E+00 
. 999603E+00 
. 999655E+00 
. 999697E+00 



200787E+00 
222538E+00 
244798E+00 
267635E+00 
291127E+00 
315344E+00 
340375E+00 
366312E+00 
393263E+00 
421337E+00 
450671E+00 
-.481411E+00 
513727E+00 
547810E+00 
583876E+00 
622189E+00 
663036E+00 
706766E+00 
753794E+00 
804590E+00 
859744E+00 
919960E+00 
986096E+00 
105921E+01 
114061E+01 
123199E+01 
133547E+01 
145388E+01 
159094E+01 
175175E+01 
194343E+01 
217624E+01 
246558E+01 
283557E+01 
332630E+01 
400966E+01 
502848E+01 
671234E+01 
100258E+02 
194057E+02 
877273E+02 
. 200011E+02 
. 102594E+02 
. 687151E+01 
.516108E+01 
. 412995E+01 
. 344001E+01 
. 294542E+01 
. 257306E+01 
. 228223E+01 
. 204848E+01 
. 185626E+01 
. 169518E+01 
. 155805E+01 
. 143975E+01 
. 133650E+01 



97 



185196E+01 

188002E+01 

190S08E+01 

193614E+01 

196420E+01 

199226E+01 

202032E+01 

204838E+01 

207644E+01 

210450E+01 

213256E+01 

216062E+01 

218868E+01 

221674E+01 

224480E+01 

227286E+01 

230092E+01 

232898E+01 

235704E+01 

238510E+01 

241316E+01 

244122E+01 

246928E+01 

24S734E+01 

252540E+01 

255346E+01 

258152E+01 

260958E+01 

263764E+01 

266570E+01 

269376E+01 

272182E+01 

274988E+01 

277794E+01 



. 999732E+00 
. 999761E+00 
. 999785E+00 
. 999805E+00 
. 999822E+00 
. 999837E+00 
. 999850E+00 
. 999861E+00 
. 999871E+00 
. 999880E+00 
. 999887E+00 
. 999894E+00 
. 999900E+00 
. 999906E+00 
. 999911E+00 
. 999915E+00 
. 999919E+00 
. 999922E+00 
. 999926E+00 
. 999929E+00 
. 999932E+00 
. 999934E+00 
. 999936E+00 
. 999938E+00 
. 999940E+00 
. 999942E+00 
. 999943E+00 
. 999945E+00 
. 999946E+00 
. 999948E+00 
. 999949E+00 
. 999950E+00 
. 999951E+00 
. 999952E+00 



. 124548E+01 
. 116453E+01 
. 109196E+01 
. 102645E+01 
. 966937E+00 
. 912553E+00 
. 862608E+00 
. 816508E+00 
. 773770E+00 
. 733993E+00 
. 696823E+00 
. 661968E+00 
. 629177E+00 
.598225E+00 
. 568932E+00 
. 541125E+00 
. 514664E+00 
.489416E+00 
.465272E+00 
. 442129E+00 
. 419903E+00 
. 398506E+00 
. 377872E+00 
. 357931E+00 
. 338631E+00 
. 319911E+00 
. 301721E+00 
. 284018E+00 
. 266771E+00 
. 249936E+00 
. 233465E+00 
. 217347E+00 
. 201533E+00 
. 186002E+00 



END OF RUN, SYSTEM # 1 
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ooooooooooonooooooooooooooonnooooonooooooooooooo 



APPENDIX B. PROGRAM AND DATA OF DIFFERENCE EQUATION 
A. PROGRAM FOR DIFFERENCE EQUATION CALCULATION 

DIFFEQ. FOR VERSION; 2/03/88 



PURPOSE: THIS PROGRAM COMPUTES THE ITERATIVE SOLUTION TO A 

LINEAR, TIME-INVARIANT (LTI) DIFFERENCE EQUATION. 

THE DIFFERENCE EQUATION MUST BE IN THE FORM: 
y(ns) = a( l)'''y(ns-l) + ... + a(N)''y(ns-N) + 

b(0)-'x(ns) + b(l)*x(ns-l) + ... + b(L)*x(ns-L). 
THE PROGRAM CONSISTS OF A MAIN PROGRAM AND T\v^0 
SUBROUTINES. SUBROUTINE diffeq IS CALLED BY THE MAIN 
PROGRAM TO ITERATIVELY SOLVE THE DIFFERENCE EQUATIONS AND 
SUBROUTINE xgen ALLOWS THE USER THE OPTION OF GENERATING 
THE INPUT SEQUENCE x() BY WRITING THE APPROPRIATE 
EQUATIONS. IF THE USER ELECTS TO GENERATE THE SEQUENCE 
x() BY USING xgen THEN THE PROGRAM MUST BE COMPILED 
AGAIN BEFORE EXECUTION. THE USER HAS THE OPTION OF 
SELECTING ONE OF TWO OPERATING MODES: BATCH OR TEST. 

IN BATCH MODE THE AMOUNT OF INTERFACE WITH THE USER 
IS MINIMIZED AND IT IS ASSUMED THAT THE INPUT DATA 
HAS BEEN STORED IN THE DEFAULT FILE 'DIFFEQ. IN'. IN 
TEST MODE THE USER IS PROMPTED FOR THE NAME OF THE 
INPUT FILE OR HAS THE OPTION OF PERFORMING A TEST RUN 
USING THE INPUT DATA STORED IN THE FILE ' DIFFEQ. TST' . 

IT IS RECOMMENDED THAT FIRST-TIME USERS SELECT THE 
TEST MODE AND MAKE A TRIAL RUN WITH THE PRE STORED 
INPUT DATA. THE TEST MODE ECHOES PORTIONS OF THE 
INPUT DATA ONTO THE MONITOR TO ALLOW VERIFICATION OF 
ITS ACCURACY. BOTH BATCH AND TEST MODES ALLOW THE 
USER TO SOLVE UP TO FOUR DIFFERENCE EQUATIONS IN A 
SINGLE PROGRAM EXECUTION. THE OUTPUT OF THE PROGRAM 
'DIFFEQ. FOR' IS STORED IN THE ARRAY y(). THE OUTPUT IS 
STORED IN TABULAR FORM IN THE OUTPUT FILE 'DIFFEQ. FOR' 

AND IN A FORM SUITABLE FOR PLOTTING IN THE FILE 
'DIFFEQ.DAT'. 












THIS PROGRAM ASSUMES THAT EACH DIFFERENCE EQUATION IS IN THE 
CANONICAL FORM; 

y(ns) = a(l)*y(ns-l) + ... + a(N)*y(ns-N) + 

b(0)*x(ns) + b(l)*x(ns-l) + ... + b(L)*x(ns-L) 

L = A NON-NEGATIVE INTEGER, THE NUMBER OF INPUT DELAYS. 

N = A NON-NEGATIVE INTEGER, THE NUMBER OF OUTPUT DELAYS. 
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C a(l)...aCN) = REAL COEFFICIENTS OF THE OUTPUT TERMS. 

C b(0)...b(L) = REAL COEFFICIENTS OF THE INPUT TERMS. 

C 

C THE INPUT PARAMETERS SHOULD BE STORED IN A FILE NAMED 
C 'DIFFEQ. IN’. ALL OF THE READ STATEMENTS USED BY THIS PROGRAM 
C REQUIRE FORMATTED INPUT. PARTICULAR ATTENTION SHOULD BE PAID 
C TO THE FORMATS, ESPECIALLY THE USE OF THE DECIMAL POINT TO 
C DENOTE ’REAL’ NUMBERS. THE INPUT PARAMETERS REQUIRED BY THE 
C PROGRAM ARE LISTED BELOW. 

C 

C 



c 

c 

c 


NAME 


TYPE 


RANGE 


(ARRAYS) 




RESTRICTIONS 


numsys 


INTEGER 






1 


<= numsys <= 4 


c 


L 


INTEGER 








0 <= L <= 128 


c 


N 


INTEGER 








0 <= N <= 128 


c 


nstop 


INTEGER 






0 


<= nstop <= 999 


c 


xsorce 


CHARACTER 








’F’ OR 's’ 


c 


b(k) 


REAL 


0,1, 


... , L 




0 <= L <= 128 


c 


a(k) 


REAL 


1,2, 


. . . , N 




0 <= N <= 128 


c 


y(k) 


REAL 


-N, . 


. . , -1 




1 <= N <= 128 


c 


x(ns) 


REAL 


0,1, .. 


. , nstop 


0 


<= nstop <= 999 



C 

C WHERE: 

C 

C numsys = THE NUMBER OF SYSTEMS TO BE EVALUATED. 

C THIS INTEGER VALUE MUST OCCUR AT THE TOP OF THE INPUT 

C FILE. IT DELINEATES THE NUMBER OF SYSTEMS TO BE READ BY 

C THE PROGRAM AND ANALYZED. FOR EACH SYSTEM (1. ..numsys) 

C THE PARAMETERS BELOW MUST APPEAR IN THE INPUT FILE. 

C 

C L = AN INTEGER VALUE THAT SPECIFIES THE MAXIMUM NUMBER OF DELAYS 
C IN THE INPUT SEQUENCE. 

C 

C N = AN INTEGER VALUE THAT SPECIFIES THE MAXIMUM NUMBER OF DELAYS 
C IN THE OUTPUT SEQUENCE. 

C 

C nstop = AN INTEGER VALUE THAT SPECIFIES THE LARGEST TIME INDEX 
C (ns) FOR WHICH THE DIFFERENCE EQUATION IS TO BE SOLVED. 

C 

C xsorce = A CHARACTER VALUE OF ’F’ OR ’S’ DENOTING WHETHER THE 
C INPUT SEQUENCE x( ) IS TO BE READ FROM THE INPUT FILE (F) 

C OR TO BE GENERATED (S) USING THE SUBROUTINE xgen. THIS 

C LATER OPTION IS ATTRACTIVE WHEN nstop IS A LARGE NUMBER 

C AND THE INPUT SEQUENCE x( ) CAN BE READILY DESCRIBED BY AN 

C ANALYTICAL EXPRESSION. IF xsorce = ’S’ THE USER MUST 

C PROVIDE THE APPROPRIATE FORTRAN STATEMENTS IN THE SPACE 

C PROVIDED IN SUBROUTINE xgen AND THE PROGRAM MUST BE 

C RECOMPILED BEFORE EXECUTION. 

C 

C b(k) = REAL COEFFICIENTS OF THE INPUT SEQUENCE x(ns-k) IN THE 
C ORDER: b(0) , b(l) , . . . , b(L). 

C 

C a(k) = REAL COEFFICIENTS OF THE OUTPUT SEQUENCE y(ns-k) IN THE 
C ORDER: a(l), a(2), ..., a(N). IF N = 0 THEN THE EQUATION 

C IS NON-RECURSIVE AND NO a(k) COEFFICIENTS SHOULD BE IN 
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THE INPUT FILE. 



C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C' 

c 

c 

c 

c 

c 



y(k) = THE INITIAL CONDITIONS FOR THE OUTPUT SEQUENCE IF THE 
DIFFERENCE EQUATION IS RECURSIVE, I.E. , N > 0. THIS 
PROGRAM CALCULATES THE SOLUTION TO THE DIFFERENCE 
EQUATION FROM ns = 0 TO ns = nstop THEREFORE THE INITIAL 
CONDITIONS y(-N) TO y(-l) MUST BE PROVIDED IN THE INPUT 
FILE IN THE ORDER: y(-N), y( -N+1) , y(-l). IF N = 0 

THEN THE EQUATION IS NON-RECURSIVE AND NO INITIAL 
CONDITIONS SHOULD BE GIVEN IN THE INPUT FILE. 

x(ns) = THE INPUT SEQUENCE. IF xsorce = 'F' THEN THE INPUT 

SEQUENCE x(0), ... , x(nstop) MUST BE PROVIDED BY THE USER 
IN THE INPUT FILE. IF xsorce = ' S' THEN THE USER HAS 
ELECTED TO GENERATE THE INPUT SEQUENCE BY PROVIDING THE 
APPROPRIATE FORTRAN STATEMENTS IN THE SUBROUTINE xgen. 

NOTE: THE INPUT FORMAT STATEMENTS OCCUR IN THE MAIN PROGRAM 

FOLLOWING THE CAPTION; ******** INPUT FORMAT ********. 
THE FORM OF THE INPUT DATA FILE IS: 



LINE # 



ENTRIES 



FORMAT 



1 

2 

NEXT NB LINES 
NEXT NA LI.NES 
NEXT NY LINES 
NEXT NX LINES 



numsys 

L,N, nstop, xsorce 
b(k) , k=0,l,. . . ,L 
a(k) , k=l,. . . ,N 
y(k) , k= -N,. . . ,-l 
x(ns), ns= 0 ,..., nstop 



il 

i3, til , i3 ,t21 , i3, t31 , al 
6fl0. 0 
6fl0. 0 
6fl0. 0 
6fl0. 0 



WHERE; NB 

NA 

NA 

NY 

NY 

NX 

NX 



1 + (L/6 ROUNDED DOWN TO THE NEXT SMALLER INTEGER) 

0 IF N = 0 OR 

1 + ((N-l)/6 ROUNDED DOWN TO THE NEXT SMALLER INTEGER) 

0 IF N = 0 OR 

1 + ((N-l)/6 ROUNDED DOWN TO THE NEXT SMALLER INTEGER) 

0 IF xsorce = 'S' OR 

1 + (nstop/6 ROUNDED DOWN TO THE NEXT SMALLER INTEGER) 
IF xsorce = 'F' 



*NOTE; FOR numsys > 1 THE FORMAT OF LINES 2. . . IS REPEATED. 



THE FORM. T f 10. 0 USED FOR INPUT DATA PERMITS THE DECIMAL 
POINT TO BE PLACED ANYWHERE IN THE FIELD OF 10 COLUMNS 
AND ALSO ALLOWS THE EXPONENTIAL FORMAT TO BE USED (EG. 
3146. 2=3. 1462E+03). 



OUTPUT ^V‘5V*‘:V‘5V^V^V^V‘5V'>V:?rVc*'5VycVr'}Wr’5V'5V*^WcVr***5V‘5V‘>V^V 



THE INPUT DATA AS WELL AS THE OUTPUT DATA ARE STORED IN TABULAR 
FORM IN THE FILE ' DIFFEQ. OUT' . ADDITIONALLY, THE INPUT SEQUENCE 
AND THE OUTPUT SEQUENCE ARE WRITTEN INTO THE FILE 'DIFFEQ.DAT' TO 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



FACILITATE PLOTTING BY A SEPARATE, USER SUPPLIED PROGRAM. THE 
FORMAT OF THE DATA IN 'DIFFEQ.DAT' IS: el2.6, 2X, el2. 6. THE FIRST 
ENTRY CORRESPONDS TO THE ORDINATE VALUE, AND THE SECOND ENTRY, THE 
ABSCISSA VALUE. ADDITIONAL HEADER INFORMATION IS WRITTEN INTO 
'DIFFEQ.DAT' TO ALLOW FOR CONTROL AND LABELING OF EACH PLOT. 



■}V**^**';V***:fWc";V***:iV**?Wc’*****V«’Vc** EXAMPLE 






THE INPUT PARAMETERS FOR THE SYSTEM DESCRIBED BELOW ARE STORED IN 
THE SAMPLE INPUT FILE 'DIFFEQ.TST' AND CAN BE USED FOR A TRIAL 
RUN IN TEST MODE. 

DIFFERENCE EQUATION: 

y(ns) = 1. 2*y(ns-l) + 1.5*x(ns) 

GOAL: TO OBTAIN THE SOLUTION TO THIS DIFFERENCE EQUATION FOR 

ns = 0 TO ns = 10, GIVEN: x(0)...x(10) = 100.0 AND 
THE INITIAL CONDITION y(-l) = 25.0. 



THE INPUT FILE IS: 



1 

000 001 

1. 5 
1 . 2 
25. 0 

100.0 100.0 

100.0 100.0 



010 F 



100.0 100.0 100.0 100.0 

100.0 100.0 100.0 



THE RESULTING OUTPUT FILE 'DIFFEQ. OUT' IS: 

INPUT DATA FOR PROBLEM Z/ 1 

PROBLEM # 1 INPUT DATA SOURCEFILE: DIFFEQ.TST 

THE NUMBER OF INPUT DELAYS: L = 0 

THE NUMBER OF OUTPUT DELAYS: N = 1 

THE VALUE OF nstop IS: 10 

THE COEFFICIENTS b(0), b(l), ..., b(L) ARE: 

. 150000E+01 



THE COEFFICIENTS a(l), ..., a(N) ARE: 
. 120000E+01 



OUTPUT DATA FOR PROBLEM # 1 

ns x(nr>) y(ns) 

-1 .OOOOcOE+00 . 250000E+02 
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0 . lOOOOOE+03 

1 . lOOOOOE+03 

2 . lOOOOOE+03 

3 . lOOOOOE+03 

4 . lOOOOOE+03 

5 . lOOOOOE+03 

6 . lOOOOOE+03 

7 . lOOOOOE+03 

8 . lOOOOOE+03 

9 . lOOOOOE+03 

10 . lOOOOOE+03 

END OF 



. 180000E+03 
. 366000E+03 
. 589200E+03 
. 857040E+03 
. 117845E+04 
. 156414E+04 
. 202697E+04 
. 258236E+04 
. 324883E+04 
. 404860E+04 
. 500832E+04 

PROBLEM # 1 



FOR ILLUSTRATIVE PURPOSES THE INPUT SEQUENCE x() COULD HAVE BEEN 
GENERATED BY SPECIFYING xsorce = 'S' AND WRITING THE APPROPRIATE 
FORTRAN STATEMENTS INTO SUBROUTINE xgen. THE STATEMENTS THAT 
COULD BE USED TO ACCOMPLISH THIS ARE WRITTEN INTO THE SUBROUTINE 
BUT ARE 'COMMENTED OUT'. 



' MA I N PROGRAM ************Vr*VrVf*Vr-VVoV*VrVrVf’;V** 



character infile*12j mode*l, xsorce*! 

real a(l: 128) , b(0:128)3 y(-128:850), x(-128:850), ii 

C PROMPT USER FOR MODE: BATCH OR TEST. 

write(* 5 1115) 
read(* 3 1117) mode 

if ((mode. eq. ’ Y’ ). or. (mode. eq. V* )) then 
mode = ^Y’ 
write(*, 1118) 
read(*, 1119) infile 
else 

infile = 'DIFFEQ. IN' 
end if 

C UNIT=1 DEFINED AS INPUT FILE. UNITS=2,3 DEFINED AS OUTPUT FILES. 

open(unit=l , f ile=inf ile, status='old' , iostat=ierr , err=999) 
open(unit=2,file=' DIFFEQ. OUT' ) 
open(unit=3 , file=' DIFFEQ. DAT' ) 

C READ INPUT PARAMETERS AND CONDUCT ERROR CHECKS. 

read( 1,1000) numsys 

if ( (numsys. le. 0) . or. (numsys. gt. 4) ) then 
write(*, 1126) numsys 

stop 'The allowed values for numsys are 1 <= numsys <= 4. ' 
end if 

numplts = 2''numsys 
writeC 3 , 2000) numplts 
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do 10 nprob=l, numsys 

read( 1,1001) L, N, nstop, xsorce 

if((L. It. 0).or. (L. gt. 128)) then 
write(* , 1124) nprob, 'L', L 

stop 'The allowed values for ' 'L' ' are: 0 <= L <= 128. ' 

endif 

if((N. It. 0). or. (N. gt. 128)) then 
write(*, 1124) nprob, 'N' , N 

stop 'The allowed values for ' 'N' ' are: 0 <= N <= 128. ' 

endif 

if((nstop. It. 0). or. (nstop. gt. 850)) then 
write(*, 1127) nprob, nstop 

stop 'The allowed values for nstop are: 0 <= nstop <= 850. ' 
endif 

if ( (xsorce. eq. ' F ' ) . or. (xsorce. eq. ' f ' ) ) then 
xsorce = 'F' 

elseif( (xsorce. eq. ' S ' ) . or. ( xsorce. eq. ' s ' ) ) then 
xsorce = 'S' 
else 

write(*, 1128) nprob, xsorce 

stop 'The allowed values for ''xsorce'' are: ''f'' or ''S''.' 

endif 

C INITIALIZE EACH ARRAY TO ZERO BEFORE EACH RUN. 

data a/128*0.0/, b/129*0. 0/ 
data y/979-’-O. 0/, x/979*0. 0/ 

C READ THE COEFFICIENTS b( ) , a( ) AND THE INITIAL CONDITIONS 
C y(-N). . .y(-l). 

read(l,1002) (b(k), k=0,L) 
if(N. gt.O) then 
icstart = -N 

readd, 1002) (a(k), k=l,N) 
read( 1,1002) (y(k), k=icstart , -1) 
endif 

C FOR xsorce = 'F' READ THE ARRAY x( ) FROM THE INPUT FILE. 

C FOR xsorce = 'S' CALL xgen TO GENERATE THE ARRAY x(). 

if(xsorce. eq. 'F' ) then 

read( 1,1002) (x(k), k=0, nstop) 
else 

call xgen(x, nstop, nprob) 
endif 

C FOR TEST MODE ECHO INPUT PARAMETERS ONTO MONITOR (UNIT = *). 
if(mode. eq. ' Y' ) then 
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write(*, 1007) 

write(*, 1120) nprob, infile 
write(*,1110) ’INPUT* 'L*. L 
write(*,1110) *0U"TUT^ *N\ N 
write(*, 1112) nstcp 
write(*, 1004) 

write(*,1005) (b(k),k=0,L) 
if(N. eq. 0) then 
write(*,1131) 
else 

write(*, 1006) 

write(*,1005) (a(k),k=l,N) 
end if 

write(*, 1123) nprob 

pause ’END OF RUN, STRIKE <CR> WHEN READY TO CONTINUE. ' 
endif 

C raiTE INPUT DATA INTO FILE: DIFFEQ. OUT. 

write( 2 , 1008) *INPUT*, nprob 
write( 2 , 1120) nprob, infile 
write(2,1110) ’INPUT*. *L*. L 
write(2,1110) ’OUTPUT^ *N\ N 
write( 2,1112) ns top 
write( 2 , 1004) 

write(2,1005) (b(k),k=0,L) 
if(N. eq. 0) then 
write(2,1131) 
else 

write( 2 , 1006) 

write(2, 1005) (a(k) ,k=l,N) 
endif 

C WRITE THE INPUT SEQUENCE INTO FILE: DIFFEQ.DAT. 

write(3,2001) nstop + 1 
write(3,^*“) ’INPUT SEQUENCE x(ns)* 
write(3U'0 ’SAMPLE // (ns)' 
write( 3 ,*) * x(ns) ’ 
do 54 ns=0, nstop 
ii = ns 

write(3,2010) ii, x(ns) 

54 continue 

C CALL diffeq TO COMPUTE THE SOLUTION TO THE DIFFERENCE EQUATION. 

call diffeq(N,L,a,b,x,y ,nstop) 

C WRITE RESULTS INTO FILE: DIFFEQ.DAT. 

write(3,2001) N + nstop + 1 
write(3,*) 'OUTPUT SEQUENCE Y(ns)' 
write(3,*) 'SAMPLE (ns)’ 
write(3,*) ’y(ns) * 
do 55 ns= -N, nstop 
ii = ns 
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write(3 ,2010) ii, y(ns) 

55 continue 

C WRITE RESULTS INTO FILE: DIFFEQ. OUT. 

write(2,I008) 'OUTPUT', nprob 
write(2 , 1129) 
do 102 ns= -N, nstop 

write(2,1130) ns, x(ns), y(ns) 

102 continue 

write( 2, 1123) nprob 

10 continue 

write(*, 1121) i 
999 close(unit=l) 
close(unit=2) 
close(unit=3) 

if ( ierr. gt. 0) then 
write(* , 1116) ierr 
end if 



C 



INPUT FORMAT ****v?*** 



1000 format(il) 

1001 format ( i3 , til , i3 , t2 1 , i3 , t31 , al) 

1002 format(6fl0. 0) 



1004 forraat(t4, 'THE COEFFICIENTS b(0), b(I), ..., b(L) ARE:',/) 

1005 format(6( lx,el2. 6)) 

1006 format(//,t4,'THE COEFFICIENTS a(l), ..., a(N) ARE:',/) 

1007 format(//////////) 

1008 format(///,tl6,a6, ' DATA FOR PROBLEM ',il,//) 

1110 format(t4, 'THE NUMBER OF ',a6,' DELAYS: ',al,' = ',i3) 

1112 format(t4, 'THE VALUE OF nstop IS: ' ,i3) 

1115 0format(lx, 'DO YOU WISH TO RUN THIS PROGRAM IN TEST', 

l' MODE ? (Y/N) <CR> : ', ,) 

1116 OformatC///, lx, 'ERROR OPENING INPUT FILE, PROGRAM TERMINATED.', 
I//, lx, 'ERROR CODE: ' ,i4, /////) 

1117 format(al) 

1118 OformatC/////, lx, 'TYPE THE NAME OF YOUR DATA FILE FOLLOWED', 

I' BY <CR>. ' ,/,' IF YOU DESIRE TO MAKE A TEST RUN USING THE' , 

2' SAMPLE DATA ALREADY STORED',/,' IN THE FILE: DIFFEQ. TST', 

3 ' TYPE: DIFFEQ. TST <CR>' ,/,' FILENAME; ', ,) 

1119 format(al2) 

1120 formate///, t4, 'PROBLEM //' ,il, ' INPUT DATA SOURCEFILE: ' ,al2) 

1121 OformatC//,' TABULAR OUTPUT DATA IS STORED IN FILE: DIFFEQ. OUT. ' 
1,/,' PLOTTING DATA IS STORED IN FILE; DIFFEQ.DAT. ') 

1122 formarCi3) 

1123 formate/, lx, 16C ’-'),’ END OF PROBLEM #', i2 ,2x, 16C '-'), //) 

1124 OformatC//,' For problem //',i2,' the value for ',al,' is: ',i3. 
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l'. This value is not allowed. ’) 

1126 formate//,' numsys = ',i4,'. This value is not allowed, ') 

1127 OformatC//,' For problem ?r',i2,' the value for ''nstop'' is: 
li3,'. This value is not allowed.') 

1128 OformatC//,' For problem ,il,' the value for ''xsorce'' is: ', 

lal,'. This value is not allow'ed. ') 

1129 format (t6, ' ns ' ,tl6 , 'x(ns) ' , t35 , ' y(ns) ' ) 

1130 format(t4,i4,tll,el4. 6,t30,el4. 6) 

1131 formate/,' THIS SYSTEM IS NON -RECURSIVE, I.E. , N = 0. ' ) 

2000 formate il) 

2001 formate i3) 

2010 formate el2. 6 , 2x, el2. 6) 

end 



C SUBROUTINE: diffeq 



C PURPOSE: THIS SUBROUTINE COMPUTES THE SOLUTION TO A DIFFERENCE 

C EQUATION. ALL PARAMETERS DESCRIBING THE EQUATION, AND 

C THE INPUT AND OUTPUT SEQUENCES xC ) AND yC) ARE PASSED 

C TO THE SUBROUTINE BY THE MAIN PROGRAM. 

subroutine diffeqCN,L, a,b,x,y ,nstop) 

real xC "128: nrtop) , yC - 128: nstop) , aCl:N), bC0:L) 

do 500 ns=0, nstop 
yCns) = 0, 0 
do 501 k=0, max(!:,L) 

yCr.s) = yCns) + aek)''yens-k) + bek)*xCns-k) 

501 continue 
500 continue 

return 

end 



C SUBROUTINE: xgen 



C PURPOSE: THIS SUBROUTINE ALLOWS THE USER TO GENERATE VALUES FOR 

C THE ARRAY xQ. IF xsorce = 'S' THE MAIN PROGRAM WILL 

C CALL THIS SUBROUTINE. IF xsorce = 'F' THIS SUBROUTINE 

C WILL NOT BE CALLED BY THE MAIN PROGRAM. 



subroutine xgenC x, nstop ,nprob) 
real xC - 128: nstop) 



C DE\T:L0P THE ALGORITHM FOR GENERATING VALUES OF xC ) IN THIS SPACE. 
C TF^ STATEMENTS TYPED IN MUST FOLLOW STANDARD FORTRAN 77 RULES AND 
C MAY USE FORTRAN 77 INTRINSIC FUNCTIONS SUCH AS; SINQ , COSQ, ... 
C NOTE THAT THE VALUE nprob CAN BE USED IN A LOGICAL 'IF' STATEMENT 
C TO MATCH THE GENERATING FUNCTIONS TO THE CORRESPONDING SYSTEM 
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EQUATION READ FROM THE INPUT FILE IF MORE THAN ONE SYSTEM OF 
EQUATIONS EXIST. AN EXAMPLE OF AN ALGORITHM GENERATING VALUES 
FOR x() IS: 



** EXAMPLE *** 

if(nprob. eq. 1) then 
do 1 k=0, nstop 
x(k) = 100. 0 
1 continue 
end if 

x(0) = 1.0 
do 1 k=l,nstop 
x(k) = 0. 0 
1 continue 






return 

end 



B. INPUT, OUTPUT DATA OF SECOND ORDER HR NOTCH FILTER: P = 0 



INPUT DATA FOR PROBLEM # 1 



PROBLEM # 1 INPUT DATA SOURCEFILE: diffeq. tst 

THE NUMBER OF INPUT DELAYS: L = 2 

THE NUMBER OF OUTPUT DELAYS: N = 2 

THE VALUE OF nstop IS: 400 

THE COEFFICIENTS b(0), b(l), ..., b(L) ARE: 

. 990000E+00 329930E+00 . 990000E+00 



THE COEFFICIENTS a(l), ..., a(N) ARE; 
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. 329930E+00 980100E+00 



OUTPUT DATA FOR PROBLEM fr 1 



ns 


x(ns) 


y(ns) 


-2 


. OOOOOOE+00 


. OOOOOOE+00 


-1 


. OOOOOOE+00 


. OOOOOOE+00 


0 


. lOOOOOE+01 


. 990000E+00 


1 


. OOOOOOE+00 


329930E-02 


2 


. OOOOOOE+00 


. 186125E-01 


3 


. OOOOOOE+00 


. 937446E-02 


4 


. OOOOOOE+00 


151492E-01 


5 


. OOOOOOE+00 


141861E-01 


6 


. OOOOOOE+00 


. 101673E-01 


7 


. OOOOOOE+00 


. 172583E-01 


8 


. OOOOOOE+00 


427095E-02 


9 


. OOOOOOE+00 


183239E-01 


10 


. OOOOOOE+00 


185966E-02 


11 


. OOOOOOE+00 


. 173457E-01 


12 


. OOOCOOE+00 


. 754553E-02 


13 


. OOOOOOE+00 


145111E-01 


14 


. OCOOOOE+00 


121830E-01 


15 


. OOOOOOE+00 


. 102028E-01 


16 


. OOOOOOE+00 


. 153068E-01 


17 


. OOOOOOE+00 


494956E-02 


18 


. OOOOOOE+00 


166352E-01 


19 


. OOOOOOE+00 


637380E-03 


20 


. OOOOOOE+00 


. 160938E-01 


21 


. OOOOOOE+00 


. 593453E-02 


22 


. OOOOOOE+00 


138156E-01 


23 


. OOCOOOE+00 


103746E-01 


24 


. OOOOOOE+00 


. 101178E-01 


25 


. OOCOOOE-rOO 


. 135063E-01 


26 


. OOCOOOE+OO 


546028E-02 


27 


. OOOOOOE+00 


150390E-01 


28 


. OOOOOOE+00 


. 389786E-03 


29 


. OOOOOOE+00 


. 14S684E-01 


30 


. OOOOOOE+00 


. 452349E-02 


31 


. OOOOOOE+00 


130801E-01 


32 


. OOOOOOE+00 


874898E-02 


33 


. OOOOOOE+00 


. 993321E-02 


34 


. OOOOOOE+00 


. 118521E-01 


35 


. OOOOOOE+00 


582516E-02 


36 


. OOOOOCE+00 


135382E-01 


37 


. OOOOOvZ+00 


. 124259E-02 


38 


. OOOOOOE+00 


. 136787E-01 


39 


. OOOOOOE+00 


. 329516E-02 


40 


. oooonoE +00 


123194E-01 


41 


. OOOOOOE+00 


729411E-02 


42 


. OOOOOOE+00 


. 966765E-02 


43 


. OOOOOOE+00 


. 103386E-01 


44 


. OOOOOOE+00 


606425E-02 


45 


. OOOOOOE+00 


121336E-01 
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46 


. OOOOOOE+00 


. 194032E-02 


47 


. OOOOOOE+00 


. 125324E-01 


48 


. OOOOOOE+00 


. 223309E-02 


49 


. OOOOOOE+00 


115462E-01 


50 


. OOOOOOE+00 


599809E-02 


51 


. OOOOOOE+00 


. 933748E-02 


52 


. OOOOOOE+00 


. 895944E-02 


53 


. OOOOOOE+00 


619567E-02 


54 


. OOOOOOE+00 


108253E-01 


55 


. OOOOOOE+00 


. 250079E-02 


56 


. OOOOOOE+00 


. 114349E-01 


57 


. OOOOOOE+00 


. 132171E-02 


58 


. OOOOOOE+00 


107713E-01 


59 


. OOOOOOE+00 


484919E-02 


60 


. OOOOOOE+00 


. 895708E-02 


61 


. OOOOOOE+00 


. 770790E-02 


62 


. OOOOOOE+00 


623577E-02 


63 


. OOOOOOE+00 


961188E-02 


64 


. OOOOOOE+00 


. 294043E-02 


65 


. OOOOOOE+00 


. 103907E-01 


66 


. OOOOOOE+00 


. 546304E-03 


67 


. OOOOOOE+00 


100037E-01 


68 


. OOOOOOE+00 


383596E-02 


69 


. OOOOOOE+00 


. 853905E-02 


70 


. OOOOOOE+00 


. 657691E-02 


71 


. OOOOOOE+00 


619920E-02 


72 


. OOOOOOE+00 


849133E-02 


73 


. OOOOOOE+00 


. 327429E-02 


74 


. OOOOOOE+00 


. 940264E-02 


75 


. OOOOOOE+00 


106918E-03 


76 


. OOOOOOE+00 


925080E-02 


77 


. OOOOOOE+00 


294733E-02 


78 


. OOOOOOE+00 


. 809430E-02 


79 


. OOOOOOE+00 


. 555923E-02 


80 


. OOOOOOE+00 


609907E-02 


81 


. OOOOOOE+00 


746086E-02 


82 


. OOOOOOE+00 


. 351613E-02 


83 


. OOOOOOE+00 


. 847247E-02 


84 


. OOOOOOE+00 


650839E-03 


85 


. OOOOOOE+00 


851860E-02 


86 


. OOOOOOE+00 


217265E-02 


87 


. OOOOOOE+00 


. 763226E-02 


88 


. OOOOOOE+00 


. 464753E-02 


89 


. OOOOOOE+00 


594701E-02 


90 


. OOOOOOE+00 


651714E-02 


91 


. OOOOOOE+00 


. 367847E-02 


92 


. OOOOOOE+00 


. 760109E-02 


93 


. OOOOOOE+00 


109744E-02 


94 


. OOOOOOE+00 


781190E-02 


95 


. OOOOOOE+00 


150178E-02 


96 


. OOOOOOE+00 


. 716096E-02 


97 


. OOOOOOE+00 


. 383451E-02 


98 


. OOOOOOE+00 


575334E-02 


99 


. OOOOOOE+00 


-.565640E-02 


100 


. OOOOOOE+00 


. 377263E-02 


101 


. OOOOOOE+00 


. 678855E-02 



no 



102 


. OOOOOOE+00 


145781E-02 


103 


. OOOOOOE+00 


713443E-02 


104 


. OOOOOOE+00 


925062E-03 


105 


. OOOOOOE+00 


. 668725E-02 


106 


. OOOOOOE+00 


. 311298E-02 


107 


. OOOOOOE+00 


552711E-02 


108 


. OOOOOOE+00 


487459E-02 


109 


. OOOOOOE+00 


. 380885E-02 


110 


. OOOOOOE+00 


. 603424E-02 


111 


. OOOOOOE+00 


174217E-02 


112 


. OOOOOOE+00 


648895E-02 


113 


. OOOOOOE+00 


433395E-03 


114 


. OOOOOOE+00 


. 621683E-02 


115 


. OOOOOOE+00 


. 247589E-02 


116 


. OOOOOOE+00 


527624E-02 


117 


. OOOOOOE+00 


416741E-02 


118 


. OOOOOOE+00 


. 379629E-02 


119 


. OOOOOOE+00 


. 533699E-02 


120 


. OOOOOOE+00 


195991E-02 


121 


. OOOOOOE+00 


587742E-02 


122 


. OOOOOOE+00 


182246E-04 


123 


. OOOOOOE+00 


. 575444E-02 


124 


. OOOOOOE+00 


. 191643E-02 


125 


. OOOOOOS+00 


500764E-02 


126 


. OOOOOOE+00 


353046E-02 


127 


. OOOOOOE+00 


. 374319E-02 


128 


. OOOOOOE+00 


. 469519E-02 


129 


. OOOOOOE+00 


211961E-02 


130 


. OOOOOOE+00 


530108E-02 


131 


. OOOOOOE+00 


. 328445E-03 


132 


. OOOOOOE+00 


. 530396E-02 


133 


. OOOOOOE+00 


. 142802E-02 


134 


. 000000E-!-00 


472726E-02 


135 


. OOOOOOE+00 


295927E-02 


136 


. OOOOOOE+00 


. 365683E-02 


137 


. OOOOOOE-rOO 


. 410688E-02 


138 


. OOOOOOE+00 


222908E-02 


139 


. OOOOOOE+00 


476059S-02 


140 


. OCOOOOE+00 


. 614057E-03 


141 


. OOOOOOE-rOO 


. 486845E-02 


142 


. OOOOOCE-l-00 


. 100441E-02 


143 


. OOOOOOE+00 


444019E-02 


144 


. OOOOOOE+00 


244937E-02 


145 


. OOOOOOE-fOO 


. 354370E-02 


146 


. OOOOOOE-i-OO 


. 356981E-02 


147 


. OOOOOOE+00 


229540E-02 


148 


. OOOOOOE-I-OO 


425609E-02 


149 


. OOOOOOE+00 


. 845509E-03 


150 


. OOOOOOE-I-OO 


. 445035E-02 


151 


. OOOOOOE-HOO 


. 639621E-03 


152 


. OOOOOOE-I-OO 


415076E-02 


153 


. OOOOOOE-I-OO 


199635E-02 


154 


. OOOOOOE-I-OO 


. 340950E-02 


155 


. OOOOOOE-I-OO 


. 308152E-02 


156 


. OOOOOOE-I-OO 


232497E-02 


157 


. OOOOOOE+00 


378727E-02 



Ill 



158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 

181 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 

201 

202 

203 

204 

205 

206 

207 

208 

209 

210 

211 

212 

213 



. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 



. 102916E-02 
.405146E-02 
. 328015E-03 
386261E-02 
159588E-02 
. 325922E-02 
. 263944E-02 
232353E-02 
335351E-02 
. 117087E-02 
. 367308E-02 
. 642924E-04 
357878E-02 
124376E-02 
. 309721E-02 
. 224087E-02 
229624E-02 
295387E-02 
. 127597E-02 
. 331607E-02 
156510E-03 
330172E-02 
935942E-03 
. 292722E-02 
. 188310E-02 
224768E-02 
258720E-02 
. 134936E-02 
. 298091E-02 
339014E-03 
303344E-02 
668554E-03 
. 275250E-02 
. 156338E-02 
218192E-02 
225215E-02 
. 139544E-02 
. 266773E-02 
487510E-03 
277549E-02 
437908E-03 
. 257578E-02 
. 127902E-02 
210253E-02 
194725E-02 
. 141823E-02 
. 237642E-02 
605957E-03 
252905E-02 
240512E-03 
. 239937E-02 
. 102735E-02 
201267E-02 
167095E-02 
. 142132E-02 
. 210663E-02 
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214 

215 

216 

217 

218 

219 

220 

221 

222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 

261 

262 

263 

264 

265 

266 

267 

268 

269 



. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOf'E+00 
. 00000C2+00 
. OOOOOOE+00 
. OOOOOOE+00 
. oonoooE+00 

. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 



697998E-03 
229500E-02 
730823E-04 
. 222522E-02 
. 805795E-03 
191508E-02 
142160E-02 
. 140794E-02 
. 185783E-02 
766969E-03 
207391E-02 
. 674610E-04 
. 205490E-02 
. 611853E-03 
181214E-02 
119756E-02 
. 138096E-02 
. 162935E-02 
815913E-03 
186612E-02 
. 183989E-03 
. 188968E-02 
. 443136E-03 
170587E-02 
997137E-03 
. 134294E-02 
. 142037E-02 
847595E-03 
167175E-02 
. 279167E-03 
. 173059E-02 
. 297362E-03 
159804E-02 
818687E-03 
. 129613E-02 
. 123003E-02 
864516E-03 
149078E-02 
. 355459E-03 
. 157839E-02 
. 172373E-03 
149011E-02 
660574E-03 
. 124251E-02 
. 105737E-02 
868928E-03 
132301E-02 
.415134E-03 
. 143365E-02 
. 661320E-04 
138330E-02 
521209E-03 
. 118381E-02 
. 901413E-03 
862852E-03 
116816E-02 
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270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 

281 

282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 



. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOO'^E+OO 
. OOOOOCE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. 0onr>Q0E+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 



. 460272E-03 
. 129677E-02 
232701E-04 
127864E-02 
399054E-03 
. 112153E-02 
. 761140E-03 
848092E-03 
102580E-02 
. 492771E-03 
. 116797E-02 
976162E-04 
117693E-02 
292632E-03 
. 105697E-02 
. 635534E-03 
826250E-03 
895491E-03 
. 514358E-03 
. 104737E-02 
158563E-03 
107885E-02 
200536E-03 
. 991213E-03 
. 523576E-03 
798744E-03 
776687E-03 
. 526597E-03 
. 934971E-03 
207643E-03 
984873E-03 
121428E-03 
. 925211E-03 
. 424267E-03 
766821E-03 
668821E-03 
. 530897E-03 
. 830670E-03 
246269E-03 
895391E-03 
540483E-04 
. 859741E-03 
. 336627E-03 
731569E-03 
571295E-03 
. 528523E-03 
. 734302E-03 
275737E-03 
810663E-03 
. 278823E-05 
. 795451E-03 
. 259710E-03 
693935E-03 
-.483492E-03 
. 520607E-03 
. 645634E-03 
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326 

327 

328 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 



. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOO-'^E+OO 
. OOOOOCE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OnncoOL+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 



297233E-03 
730852E-03 
. 501878E-04 
. 732867E-03 
. 192606E-03 
654736E-03 
404790E-03 
. 508155E-03 
.564390E-03 
311833E-03 
656042E-03 
. 891797E-04 
. 672410E-03 
. 134443E-03 
614672E-03 
334566E-03 
. 492056E-03 
. 490253E-03 
320515E-03 
586244E-03 
. 120718E-03 
. 614406E-03 
. 843958E-04 
574335E-03 
272207E-03 
. 473097E-03 
. 422879E-03 
324162E-03 
521414E-03 
. 145681E-03 
. 559102E-03 
. 416829E-04 
534224E-03 
217110E-03 
. 451961E-03 
. 361905E-03 
323564E-03 
461457E-03 
. 164877E-03 
. 506671E-03 
. 557033E-05 
494751E-03 
168693E-03 
.429248E-03 
. 306958E-03 
319432E-03 
406239E-03 
. 179045E-03 
. 457227E-03 
246287E-04 
456254E-03 
126393E-03 
•405474E-03 
. 257656E-03 
312396E-03 
355598E-03 
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382 


. OOOOOOE+00 


383 


. OOOOOOE+00 


384 


. OOOOOOE+00 


385 


. OOOOOOE+00 


386 


. OOOOOOE+00 


387 


. OOOOOOE+00 


388 


. OOOOOOE+00 


389 


. OOOOOOE+00 


390 


. OOOOOOE+00 


391 


. OOOOOOE+00 


392 


. OOOOOOE+00 


393 


. OOOOOOE+00 


394 


. OOOOOOE+00 


395 


. OOOOOOE+00 


396 


. OOOOOOE+00 


397 


. OOOOOOE+00 


398 


. OOOOOOE+00 


399 


. OOOOOOE+00 


400 


. OOOOOOE+00 




END OF 



. 188857E-03 
.410831E-03 
495535E-04 
419005E-03 
896748E-04 
. 381080E-03 
. 213620E-03 
303017E-03 
309343E-03 
. 194925E-03 
. 367499E-03 
697972E-04 
383214E-03 
580256E-04 
. 356444E-03 
. 174472E-03 
291787E-03 
267270E-03 
. 197800E-03 

PROBLEM //I 



C. INPUT, OUTPUT DATA OF THE AUGMENTED HR NOTCH FILTER: P= 1 



INPUT DATA FOR PROBLEM 1 



PROBLEM if: 1 INIUT DATA SOURCEFILE: diffeq. tst 

THE NUMBER OF INPUT DELAYS: L = 3 

THE NUMBER OF OUTPUT DELAYS: N = 3 

THE VALUE OF nstop IS: 400 

THE COEFFICIENTS b(0), b(l), b(L) ARE: 

. 990000E+00 329900E-02 .881146E+00 . 326631E+00 



THE COEFFICIENTS a(l), a(N) ARE: 

. OOOOOOE+00 -.871246E+00 323364E+00 



OUTPUT DATA FOR PROBLEM # 1 
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ns 

-3 

-2 

-1 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 



x(ns) 

. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. lOOOOOE+01 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOCE+00 
. OOOOOCE+00 
. OOOOOOE+OO 
. OOOOOOE+OO 
. ooonooE+00 

. OOOOOOE+OO 
. OOOOOOE+OO 



y(ns) 

. OOOOOOE+OO 
. OOOOOOE+OO 
. OOOOOOE+OO 
. 990000E+00 
329900E-02 
. 186127E-01 
, 937460E-02 
151494E-01 
141862E-01 
. 101675E-01 
. 172585E-01 
427105E-02 
183242E-01 
185964E-02 
. 173460E-01 
. 754559E-02 
145113E-01 
121831E-01 
. 102029E-01 
. 153069E-01 
494966E-02 
166353E-01 
637331E-03 
. 160940E-01 
. 593455E-02 
138158E-01 
103747E-01 
. 101179E-01 
. 135064E-01 
546039E-02 
150392E-01 
. 389857E-03 
. 148685E-01 
. 452347E-02 
130802E-01 
874900E-02 
. 993335E-02 
. 118522E-01 
582528E-02 
135383E-01 
. 124268E-02 
. 136788E-01 
. 329511E-02 
123195E-01 
729410E-02 
. 966777E-02 
. 103386E-01 
606436E-02 
121337E-01 
. 194041E-02 
. 125324E-01 
. 223302E-02 
115463E-01 
599805E-02 
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51 


. OOOOOOE+00 


. 933758E-02 


52 


. OOOOOOE+00 


. 895943E-02 


53 


. OOOOOOE+00 


619577E-02 


54 


. OOOOOOE+00 


108253E-01 


55 


. OOOOOOE+00 


. 250088E-02 


56 


. OOOOOOE+00 


. 114350E-01 


57 


. OOOOOOE+00 


. 132163E-02 


58 


. OOOOOOE+00 


107714E-01 


59 


. OOOOOOE+00 


484913E-02 


60 


. OOOOOOE+00 


. 895716E-02 


61 


. OOOOOOE+00 


. 770786E-02 


62 


. OOOOOOE+00 


623586E-02 


63 


. OOOOOOE+00 


961187E-02 


64 


. OOOOOOE+00 


. 294052E-02 


65 


. OOOOOOE+00 


. 103908E-01 


66 


. OOOOOOE+00 


. 546215E-03 


67 


. OOOOOOE+00 


100038E-01 


68 


. OOOOOOE+00 


383588E-02 


69 


. OOOOOOE+00 


. 853911E-02 


70 


. OOOOOOE+00 


. 657685E-02 


71 


. OOOOOOE+00 


619928E-02 


72 


. OOOOOOE+00 


849130E-02 


73 


. OOOOOOE+00 


. 327438E-02 


74 


. OOOOOOE+00 


. 940263E-02 


75 


. OOOOOOE+00 


107009E-03 


76 


. OOOOOOE+00 


925082E-02 


77 


. OOOOOOE+00 


294724E-02 


78 


. OOOOOOE+00 


. 809434E-02 


79 


. OOOOOOE+00 


. 555915E-02 


80 


. OOOOOOE+00 


609913E-02 


81 


. OOOOOOE+00 


746081E-02 


82 


. OOOOOOE+00 


. 351621E-02 


83 


. OOOOOOE+00 


. 847244E-02 


84 


. OOOOOOE+00 


650931E-03 


85 


. OOOOOOE+00 


851860E-02 


86 


. OOOOOOE+00 


217256E-02 


87 


. OOOOOOE+00 


. 763228E-02 


88 


. OOOOOOE+00 


. 464744E-02 


89 


. OOOOOOE+00 


594707E-02 


90 


. OOOCOOE+00 


651707E-02 


91 


. OOOOOOE+00 


. 367854E-02 


92 


. OOOOOOE+00 


. 760104E-02 


93 


. OOOOOOE+00 


109753E-02 


94 


. OOOOOOE+00 


781188E-02 


95 


. OOOOOOE+00 


150168E-02 


96 


. OOOOOOE+00 


. 716097E-02 


97 


. OOOOOOE+00 


. 383442E-02 


98 


. OOOOOOE+00 


575338E-02 


99 


. OOOOOOE+00 


565632E-02 


100 


. OOOOOOE+00 


. 377269E-02 


101 


. OOOOOOE+00 


, 678848E-02 


102 


. OOOOOOE+00 


145789E-02 


103 


. OOOOOOE+00 


713439E-02 


104 


. OOOOOOE+00 


924967E-03 


105 


. OOOOOOE+00 


. 668724E-02 


106 


. OOOOOOE+00 


. 311288E-02 
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107 

108 

109 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 



. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+CO 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 



552713E-02 
-.487450E-02 
. 380890E-02 
. 603416E-02 
174225E-02 
648890E-02 
-.433302E-03 
. 621681E-02 
. 247579E-02 
527625E-02 
-.416731E-02 
. 379633E-02 
. 533691E-02 
195998E-02 
587736E-02 
181364E-04 
. 575441E-02 
. 191633E-02 
500764E-02 
353036E-02 
. 374322E-02 
. 469510E-02 
211967E-02 
530101E-02 
. 328528E-03 
. 530391E-02 
. 142793E-02 
472725E-02 
295917E-02 
. 365685E-02 
. 410679E-02 
222913E-02 
476052E-02 
. 614133E-03 
. 486840E-02 
. 100432E-02 
444016E-02 
244927E-02 
. 354371E-02 
. 356971E-02 
229544E-02 
-.425600E-02 
. 845577E-03 
.445029E-02 
. 639532E-03 
415072E-02 
199625E-02 
. 340950E-02 
. 308142E-02 
232500E-02 
378719E-02 
. 102922E-02 
. 405139E-02 
. 327932E-03 
386257E-02 
159578E-02 
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163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 

181 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 

201 

202 

203 

204 

205 

206 

207 

208 

209 

210 

211 

212 

213 

214 

215 

216 

217 

218 



. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOCE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 



. 325921E-02 
. 263934E-02 
232355E-02 
335342E-02 
. 117092E-02 
. 367301E-02 
. 642165E-04 
357873E-02 
124367E-02 
. 309719E-02 
. 224077E-02 
229626E-02 
295378E-02 
. 127602E-02 
. 331600E-02 
156579E-03 
330167E-02 
935856E-03 
. 292720E-02 
. 188300E-02 
224769E-02 
258711E-02 
. 134939E-02 
. 298083E-02 
339076E-03 
303338E-02 
668474E-03 
. 275246E-02 
. 156329E-02 
218191E-02 
225206E-02 
. 139547E-02 
. 266765E-02 
487564E-03 
277542E-02 
437834E-03 
. 257574E-02 
. 127893E-02 
210252E-02 
194717E-02 
. 141825E-02 
. 237634E-02 
606003E-03 
252899E-02 
240445E-03 
. 239933E-02 
. 102727E-02 
201266E-02 
167086E-02 
. 142134E-02 
. 210655E-02 
698037E-03 
229493E-02 
730215E-04 
. 222517E-02 
. 805719E-03 
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219 

220 

221 

222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 

261 

262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 



. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOCE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOCE+00 
. OOOOOLE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOnooOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOCE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOCE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOCOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 



191506E-02 
142152E-02 
. 140795E-02 
. 185775E-02 
767001E-03 
207384E-02 
. 675152E-04 
. 205485E-02 
. 611783E-03 
181211E-02 
119748E-02 
. 138096E-02 
. 162927E-02 
815938E-03 
186605E-02 
. 184036E-03 
. 188963E-02 
. 443071E-03 
170584E-02 
997063E-03 
. 134294E-02 
. 142030E-02 
847613E-03 
167168E-02 
. 279207E-03 
. 173054E-02 
. 297304E-03 
159801E-02 
818618E-03 
. 129612E-02 
. 122996E-02 
864528E-03 
149071E-02 
. 355493E-03 
. 157833E-02 
. 172321E-03 
149007E-02 
660510E-03 
. 124250E-02 
. 105730E-02 
868935E-03 
132295E-02 
. 415163E-03 
. 143360E-02 
. 660854E-04 
138326E-02 
521150E-03 
. 118379E-02 
. 901347E-03 
862854E-03 
116809E-02 
. 460295E-03 
. 129671E-02 
233112E-04 
127860E-02 
399000E-03 
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275 


. OOOOOOE+00 


. 112151E-02 


276 


. OOOOOOE+00 


. 761079E-03 


111 


. OOOOOOE+00 


848089E-03 


278 


. OOOOOOE+00 


102574E-02 


279 


. OOOOOOE+00 


. 492789E-03 


280 


. OOOOOOE+00 


. 116792E-02 


281 


. OOOOOOE+00 


976517E-04 


282 


. OOOOOOE+00 


117689E-02 


283 


. OOOOOOE+00 


292583E-03 


284 


. OOOOOOE+00 


. 105694E-02 


285 


. OOOOOOE+00 


. 635477E-03 


286 


. OOOOOOE+00 


826244E-03 


287 


. OOOOOOE+00 


895433E-03 


288 


. OOOOOOE+00 


. 514371E-03 


289 


. OOOOOOE+00 


. 104732E-02 


290 


. OOOOOOE+00 


158593E-03 


291 


. OOOOOOE+00 


107880E-02 


292 


. OOOOOOE+00 


200492E-03 


293 


. OOOOOOE+00 


. 991185E-03 


294 


. OOOOOOE+00 


. 523523E-03 


295 


. 00000’'E+00 


798734E-03 


296 


. OOOOOCE+00 


776631E-03 


297 


. OOOOOOE+00 


. 526606E-03 


298 


. OOOOOOE+00 


. 934919E-03 


299 


. OOnocOE+00 


207668E-03 


300 


. OOOOOOE+00 


984830E-03 


301 


. OOOOOOE+00 


121389E-03 


302 


. OOOOOOE+00 


. 925181E-03 


303 


. OOOOOOE+00 


. 424218E-03 


304 


. OOOOOOE+00 


766808E-03 


305 


. OOOOOOE+00 


668769E-03 


306 


. OOOOOOE+00 


. 530901E-03 


307 


. OOOOOOE+00 


. 830620E-03 


308 


. OOOOOOE+00 


246290E-03 


309 


. OOOOOOE+00 


895349E-03 


310 


. OOOOOOE+00 


540135E-04 


311 


. OOOOOOE+00 


. 859710E-03 


312 


. OOOOOOE+00 


. 336583E-03 


313 


. OOOOOOE+00 


731553E-03 


314 


. OOOOOOE+00 


571245E-03 


315 


. OOOOOOE+00 


. 528524E-03 


316 


. OOOOOOE+00 


. 734253E-03 


317 


. OOOOOOE+00 


275754E-03 


318 


. OOOOOOE+00 


810621E-03 


319 


. OOOOOOE+00 


. 281864E-05 


320 


. OOOOOOE+00 


. 795419E-03 


321 


. OOOOOOE+00 


. 259670E-03 


322 


. OOOOOOE+00 


693917E-03 


323 


. OOOOOOE+00 


483446E-03 


324 


. OOOOOOE+00 


. 520605E-03 


325 


. OOOOOOE+00 


. 645588E-03 


326 


. OOOOOOE+00 


297246E-03 


327 


. OOOOOOE+00 


730811E-03 


328 


. OOOOOOE+00 


. 502140E-04 


329 


. OOOOOOE+00 


. 732835E-03 


330 


. OOOOOOE+00 


. 192569E-03 
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331 

332 

333 

334 

335 

336 

337 

338 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 



. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOCE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOCOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOO' E+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOF.+OO 
. OCCCOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOCE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 



654717E-03 
404747E-03 
. 508149E-03 
. 564346E-03 
311842E-03 
656002E-03 
. 892021E-04 
. 672377E-03 
. 134410E-03 
614651E-03 
334527E-03 
.492049E-03 
. 490211E-03 
320521E-03 
586206E-03 
. 120736E-03 
. 614374E-03 
. 843669E-04 
574313E-03 
272171E-03 
. 473087E-03 
. 422840E-03 
324164E-03 
521377E-03 
. 145696E-03 
. 559070E-03 
. 416575E-04 
534201E-03 
217077E-03 
. 451950E-03 
. 361869E-03 
323564E-03 
461421E-03 
. 164889E-03 
. 506640E-03 
. 554827E-05 
494727E-03 
168663E-03 
. 429235E-03 
. 306924E-03 
319430E-03 
406206E-03 
. 179054E-03 
. 457197E-03 
246476E-04 
456231E-03 
126367E-03 
. 405459E-03 
. 257625E-03 
312392E-03 
355566E-03 
. 188864E-03 
. 410802E-03 
495696E-04 
418981E-03 
896512E-04 
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387 


. OOOOOOE+00 


388 


. OOOOOOE+00 


389 


. OOOOOOE+00 


390 


. OOOOOCE+00 


391 


. OOOOOOE+00 


392 


. OOOOOOE+00 


393 


. OOOOOOE+00 


394 


. OOOOOOE+00 


395 


. OOOOOOE+00 


396 


. OOOOOOE+00 


397 


. OOOOOOE+00 


398 


. OOOOOOE+00 


399 


. OOOOOOE+00 


400 


. OOOOOOE+00 




END OF 



. 381065E-03 
. 213592E-03 
303011E-03 
309313E-03 
. 194929E-03 
. 367471E-03 
698106E-04 
383191E-03 
580047E-04 
. 356428E-03 
. 174446E-03 
291779E-03 
267242E-03 
. 197802E-03 

PROBLEM #1 



D. INPUT, OUTPUT DATA OF THE AUGMENTED HR NOTCH FILTER: P 



INPUT DATA FOR PROBLEM // 1 



PROBLEM // 1 INPUT DATA SOURCEFILE: dif feq. tst 

THE NUMBER OF INPUT DELAYS: L = 4 

THE NUMBER OF OUTPUT DELAYS: N = 4 

THE VALUE OF nstop IS: 400 

THE COEFFICIENTS b(0), b(l), b(L) ARE: 

. 990000E+00 -.329900E-02 . 186124E-01 . 614081E+00 862533E+00 



THE COEFFICIENTS a(l), a(N) ARE: 

. OOOOOOE+00 . nnooOOE+00 610814E+00 . 853908E+00 



OUTPUT DATA FOR PROBLEM // 1 



ns 


x(ns) 


y(ns) 


-4 


. OOOOOOE+00 


. OOOOOOE+00 


-3 


. OOOOOOE+00 


. OOOOOOE+00 


-2 


. OOOOOOE+00 


. OOOOOOE+00 
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-1 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 



. OOOOOOE+00 
. lOOOOOE+01 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOf'^E+00 
. OOOOOCE+00 
. OOOOOOE+00 
. OOOOOOE+00 



. OOOOOOE+00 
. 990000E+00 
329900E-02 
. 186124E-01 
. 937506E-02 
151490E-01 
141858E-01 
. 101669E-01 
. 172586E-01 
427097E-02 
183234E-01 
186026E-02 
. 173461E-01 
. 754517E-02 
145102E-01 
121837E-01 
. 102032E-01 
. 153059E-01 
494841E-02 
166360E-01 
636442E-03 
. 160924E-01 
. 593604E-02 
138169E-01 
103729E-01 
. 101156E-01 
. 135084E-01 
546244E-02 
150363E-01 
. 386692E-03 
. 148715E-01 
. 451996E-02 
130758E-01 
875350E-02 
. 993800E-02 
. 118465E-01 
581878E-02 
135450E-01 
. 125013E-02 
. 136700E-01 
. 330475E-02 
123297E-01 
728235E-02 
. 965436E-02 
. 103531E-01 
608030E-02 
121155E-01 
. 192010E-02 
. 125545E-01 
. 220828E-02 
115183E-01 
602890E-02 
. 937158E-02 
. 892123E-02 
615305E-02 
108724E-01 
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55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 



. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00. 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 



. 255325E-02 
. 113763E-01 
. 138688E-02 
108436E-01 
476854E-02 
. 886717E-02 
. 780770E-02 
634675E-02 
948809E-02 
. 280269E-02 
. 105437E-01 
. 375915E-03 
981388E-02 
404702E-02 
. 877377E-02 
. 631545E-02 
590817E-02 
881493E-02 
. 363443E-02 
. 900161E-02 
. 339247E-03 
974710E-02 
239484E-02 
. 747933E-02 
. 624335E-02 
686032E-02 
661345E-02 
. 257313E-02 
. 952163E-02 
181849E-02 
721898E-02 
361872E-02 
. 924135E-02 
. 285663E-02 
395398E-02 
873480E-02 
. 614640E-02 
. 485444E-02 
. 195900E-02 
112130E-01 
. 228329E-02 
. 294866E-02 
. 852188E-02 
109696E-01 
. 148640E-03 
268740E-02 
. 139773E-01 
945779E-02 
. 176843E-02 
108323E-01 
. 177122E-01 
915626E-02 
. 812659E-02 
200687E-01 
. 207174E-01 
127824E-01 
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Ill 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

1A4 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 



. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 



. 191976E-01 
297913E-01 
. 254984E-01 
226412E-01 
. 345899E-01 
410138E-01 
. 356029E-01 
404615E-01 
. 545884E-01 
567688E-01 
. 551160E-01 
678938E-01 
. 812886E-01 
821409E-01 
. 885345E-01 
107627E+00 
. 119586E+00 
124219E+00 
. 141341E+00 
164948E+00 
. 177990E+00 
192404E+00 
. 221445E+00 
249570E+00 
. 269510E+00 
299557E+00 
. 341534E+00 
377730E+00 
. 413111E+00 
464408E+00 
. 522361E+00 
574880E+00 
. 636425E+00 
715627E+00 
. 797194E+00 
879632E+00 
. 980564E+00 
109802E+01 
. 121802E+01 
135007E+01 
. 150800E+01 
168159E+01 
. 186472E+01 
207394E+01 
. 231483E+01 
257492E+01 
. 285909E+01 
318488E+01 
. 354945E+01 
394512E+01 
. 438677E+01 
488765E+01 
. 544063E+01 
604826E+01 
. 673134E+01 
749681E+01 
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167 


. OOOOOOE+00 


. 834016E+01 


168 


. OOOOO'-E+OO 


927626E+01 


169 


. OOOOOCE+00 


. 103271E+02 


170 


. OOOOOOE+00 


114959E+02 


171 


. OOOOOOE+00 


. 127878E+02 


172 


. OOOOCOE+00 


142290E+02 


173 


. OOOOOOE+00 


. 158402E+02 


174 


. OOOOOOE+00 


176274E+02 


175 


. OOOOOOE+00 


. 196109E+02 


176 


. OOOOOOE+00 


218257E+02 


111 


. OOOOOOE+00 


. 242932E+02 


178 


. OOOOOOE+00 


270308E+02 


179 


. OOOOOOE+00 


. 300773E+02 


180 


. OOOOOOE+00 


334757E+02 


181 


. OOOOOOE+00 


. 372549E+02 


182 


. OOOOOOE+00 


414534E+02 


183 


. OOOOOOE+00 


. 461307E+02 


184 


. OOOOOOE+00 


513410E+02 


185 


. OOOOOOE+00 


. 571326E+02 


186 


. OOOOOOE+00 


635747E+02 


187 


. OOOOOOE+00 


. 707512E+02 


188 


. OOOOOOE+00 


787379E+02 


189 


. OOOOOOE+00 


. 876183E+02 


190 


. OOOOOOE+00 


975028E+02 


191 


. OOOOOOE+00 


. 108509E+03 


192 


. OOOOOOE+00 


120753E+03 


193 


. OOOOOOE+00 


. 134374E+03 


194 


. OOOOOOE+00 


149537E+03 


195 


. OOOOOOE+00 


. 166415E+03 


196 


. OOOOOOE+00 


185190E+03 


197 


. OOOOOOE+00 


. 206083E+03 


198 


. OOOOOOE+00 


229340E+03 


199 


. OOOOOOE+00 


. 255219E+03 


200 


. OOOOOOE+00 


284013E+03 


201 


. OOOOOOE+00 


. 316059E+03 


202 


. OOOOOOE+00 


351727E+03 


203 


. OOOOOOE+00 


. 391413E+03 


204 


. OOOOOOE+00 


435575E+03 


205 


. OOOOOOE+00 


. 484725E+03 


206 


. OOOOOOE+00 


539423E+03 


207 


. OOOOOOE+00 


. 600286E+03 


208 


. OOOOOOE+00 


668018E+03 


209 


. OOOOOOE+00 


. 743398E+03 


210 


. OOOOOOE+00 


827281E+03 


211 


. OOOOOOE+00 


. 920624E+03 


212 


. OOOOOOE+00 


102450E+04 


213 


. OOOOOOE+00 


. 114011E+04 


214 


. OOOOOOE+00 


126875E+04 


215 


. OOOOOOE+00 


. 141191E+04 


216 


. OOOOOOE+00 


157123E+04 


217 


. OOOOOOE+00 


. 174852E+04 


218 


. OOOOOOE+00 


194581E+04 


219 


. OOOOOOE+00 


. 216537E+04 


220 


. OOOOOOE+00 


240970E+04 


221 


. OOOOOOE+00 


. 268160E+04 


222 


. OOOOOOE+00 


298418E+04 
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223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 

261 

262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 



. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOCE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OnoCOGL+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 



. 332090E+04 
369562E+04 
. 411262E+04 
457667E+04 
. 509308E+04 
566777E+04 
. 630729E+04 
701898E+04 
. 781098E+04 
869233E+04 
. 967314E+04 
107646E+05 
. 119793E+05 
133309E+05 
. 14835 lE+05 
165091E+05 
. 183719E+05 
204449E+05 
. 227518E+05 
253191E+05 
. 281760E+05 
313552E+05 
. 348932E+05 
388304E+05 
. 432119E+05 
480877E+05 
. 535138E+05 
595520E+05 
. 662716E+05 
737494E+05 
. 820710E+05 
913316E+05 
. 101637E+06 
113105E+06 
. 125868E+06 
140070E+06 
. 155875E+06 
173463E+06 
. 193036E+06 
214818E+06 
. 239057E+06 
266031E+06 
. 296049E+06 
329454E+06 
. 366628E+06 
407997E+06 
. 454033E+06 
505265E+06 
. 562277E+06 
625722E+06 
, 696325E+06 
774896E+06 
. 862332E+06 
959634E+06 
. 106792E+07 
11834 lE+07 
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279 

280 

281 

282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 

334 



. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
, OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 

. ooooo;'E+oo 

. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOCGOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 



. 132251E+07 
147174E+07 
. 163780E+07 
182260E+07 
. 202826E+07 
225712E+07 
. 251180E+07 
279523E+07 
. 311063E+07 
346162E+07 
. 385221E+07 
-.428688E+07 
. 477059E+07 
530889E+07 
. 590792E+07 
657455E+07 
. 731639E+07 
814194E+07 
. 906065E+07 
100830E+08 
. 112207E+08 
124868E+08 
. 138958E+08 
154637E+08 
. 172086E+08 
191504E+08 
. 213112E+08 
237159E+08 
. 263919E+08 
293698E+08 
. 326838E+08 
363717E+08 
. 404757E+08 
450429E+08 
. 501253E+08 
557812E+08 
. 620754E+08 
690797E+08 
. 768744E+08 
855485E+08 
. 952015E+08 
105944E+09 
. 117898E+09 
131201E+09 
. 146005E+09 
162480E+09 
. 180813E+09 
201216E+09 
. 223920E+09 
249186E+09 
. 277303E+09 
308593E+09 
. 343413E+09 
382163E+09 
. 425284E+09 
473272E+09 
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335 

336 

337 

338 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

35 1 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 



. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOCOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOO^'E+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE-^00 
. OCOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 



. 526674E+09 
586101E+09 
. 652235E+09 
725830E+09 
. 807730E+09 
898871E+09 
. 100030E+10 
111316E+10 
. 123877E+10 
137855E+10 
. 153410E+10 
170720E+10 
. 189983E+10 
211420E+10 
. 235276E+10 
261823E+10 
. 291366E+10 
324243E+10 
. 360829E+10 
401544E+10 
.446852E+10 
-.497273E+10 
. 553384E+10 
615825E+10 
. 685312E+10 
762640E+10 
. 848693E+10 
944456E+10 
. 105102E+11 
116962E+11 
. 130159E+11 
144846E+11 
. 161190E+11 
179378E+11 
. 199618E+11 
222142E+11 
. 247208E+11 
275102E+11 
. 306143E+11 
340687E+11 
. 379129E+11 
421908E+11 
. 469514E+11 
522492E+11 
. 581448E+11 
647056E+11 
. 720067E+11 
801317E+11 
. 891734E+11 
992354E+11 
. 110433E+12 
122893E+12 
. 136760E+12 
152192E+12 
. 169364E+12 
188475E+12 
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391 


. OOOOOOE+00 


. 209741E+12 


392 


. OOOOOOE+00 


233408E+12 


393 


. OOOOOOE+00 


. 259745E+12 


394 


. OOOOOOE+00 


289053E+12 


395 


. OOOOOOE+00 


. 321669E+12 


396 


. OOOOOOE+00 


357964E+12 


397 


. OOOOOOE+00 


. 398356E+12 


398 


. OOOOOOE+00 


443305E+12 


399 


. OOOOOOE+00 


.493325E+12 


400 


. OOOOOOE+OO 


548990E+12 




END OF 


PROBLEM #1 



E. INPUT, OUTPUT DATA OF THE AUGMENTED HR NOTCH FILTER: P = 3 



INPUT DATA FOR PROBLEM 1 



PROBLEM # 1 INPUT DATA SOURCEFILE: dif feq. tst 

THE NUMBER OF INPUT DELAYS: L = 5 

THE NUMBER OF OUTPUT DELAYS: N = 5 

THE VALUE OF nstop IS; 400 

THE COEFFICIENTS b(0), b( 1) , b(L) ARE: 



990000E+00 


329900E-02 


. 186124E-01 


. 937440E-02 


-.661007E+00 -.604706E+00 


THE COEFFICIENTS a(l), . 


. . , a(N) ARE: 






OOOOOOE+OO 


. OOOOOOE+OO 


. OOOOOOE+OO 


. 652382E+00 


. 598659E+00 



OUTPUT DATA FOR PROBLEM # 1 



ns 


xCns) 


y(ns) 


-5 


. OOOOOOE+OO 


. OOOOOOE+OO 


-4 


. OOOOOOE+OO 


. OOOOOOE+OO 


-3 


. OOOOOOE+OO 


. OOOOOOE+OO 


-2 


. OOOOOOE+OO 


. OOOOOOE+OO 


-1 


. OOOOOOE+OO 


. OOOOOOE+OO 


0 


. lOOOOOE+01 


. 990000E+00 


1 


. OOOOOOE+OO 


329900E-02 
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2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 



. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 



. 186124E-01 
. 937440E-02 
151485E-01 
141854E-01 
. 101674E-01 
. 172582E-01 
427053E-02 
183231E-01 
185917E-02 
. 173458E-01 
. 754575E-02 
145103E-01 
121822E-01 
. 102031E-01 
. 153069E-01 
-.494890E-02 
16634 lE-01 
636683E-03 
. 160941E-01 
. 593505E-02 
138145E-01 
103735E-01 
. 101184E-01 
. 135068E-01 
-.545927E-02 
150377E-01 
. 390820E-03 
. 148691E-01 
.452445E-02 
130786E-01 
874750E-02 
. 993428E-02 
. 118532E-01 
582363E-02 
135363E-01 
. 124418E-02 
. 136801E-01 
. 329678E-02 
123172E-01 
729197E-02 
. 966946E-02 
. 103405E-01 
606189E-02 
121310E-01 
. 194278E-02 
. 125346E-01 
. 223574E-02 
115430E-01 
599489E-02 
. 934044E-02 
. 896254E-02 
619203E-02 
108213E-01 
. 250464E-02 
. 114387E-01 
. 132593E-02 
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58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 



. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
. OOOOOOE+00 
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